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^k'*The  probability  of  survival  Is  predicted  of  people  located  In  conventional » 
expediently  upgraded  basements  when  subjected  to  the  blast  effects  produced  by 
the  detonation  of  a  1-MT  weapon  near  the  ground  surface.  Two  categories  of 
potential  shelters  are  considered  here,  1.e.f  engineered  buildings  and  single* 
family  residences. 

The  first  category  Included  12  basements  designed  for  live  loads  In  the 
range  from  50  psf  to  250  psf  and  slab  spans  from  12  ft  to  20  ft.  Each  of  these 
was  analyzed  as  expediently  upgraded  using  four  different  upgrading  schemes. 

An  expedient  upgrading  scheme  involves  strengthening  the  slab  over  the  basement 
by  providing  Intermediate  supports  and  blocking  off  all  openings  Into  the  base* 
ment.  This  resulted  In  60  shelters  of  different  strengths  which  Include  the 
conventional,  unupgraded  slabs  as  base  cases. 

The  second  category  Included  four  conventional  single-family  dwellings  with 
full  basements.  Each  was  evaluated  when  upgraded  using  at"studwali,,,uMrad1ng 
concept.  Two  of  the  basements  were  reevaluated  using  the  "post  and  beam*' up¬ 
grading  concept.  These  upgrading  concepts  are  essentially  similar  and  were  used 
as  Intermediate  supports  for  strengthening  the  joist  floor  systems. 

A  probability  of  survival  function  was  developed  far  each  shelter  and  each 
upgrading  scheme.  “The  procedure  used  to  accomplish  this  consists  of  two  parts. 
The  first  Is  a  probabilistic  structural  analysis  which  determines  the  probability 
of  failure  to  the  shelter  envelope.  The  second  Is  a  probabilistic  people  sur¬ 
vival  analysis -which  considers  two  casualty  producing  mechanisms,  l.e.,  debris 
effects  fresTthe  collapse  of  the  overhead  slab  and  primary  blast.  The  probabil¬ 
ity  of  structural  failure  Is  made  use  of  In  computing  the  probability  of  survival 
against  debris  effects. 

>The  report  Includes  a  description  of  the  shelters  analyzed,  a  description  of 
the  method  used  in  performing  the  analysis,  detailed  results,  conclusions 
and  recommendations. 
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1 .  INTRODUCTION 


Eventually,  existing  buildings  may  need  to  be  surveyed  and  designated  for 
upgrading  to  provide  fallout  and  blast  protection  for  evacuees  and  workers  of 
critical  Industries  In  the  event  of  a  nuclear  weapon  attack. 

This  study  was  concerned  with  predicting  the  probability  of  survival  of 
people  located  In  expediently  upgraded  conventional  basements  when  subjected 
to  the  blast  effects  of  a  1-MT  weapon  detonated  near  the  ground  surface.  Two 
categories  of  basements  are  considered,  l.e.,  basements  of  engineered  buildings 
and  basements  of  single  family  residences. 

The  first  category  refers  to  low-rise  engineered  buildings  with  basements. 
The  first  floor  Is  at  grade  and,  therefore,  the  slab  over  the  basement  Is 
directly  exposed  to  the  blast  load.  The  basement  walls,  by  virtue  of  their 
location,  are  not  directly  exposed  to  the  blast.  In  the  analysis  performed, 
the  first  floor  slab  Is  treated  as  the  primary  structural  component.  Its 
collapse  will  result  In  casualties  due  to  debris  Impact  and  due  to  blast  loads 
entering  shelter  areas  when  the  shelter  envelope  Is  breached.  Interior  and 
exterior  basement  walls  are  assumed  to  be  stronger  than  the  overhead  slab  and 
are  not  explicitly  considered  In  the  analysis. 

The  first  floor  slab  can  be  most  expediently  strengthened  by  reducing  Its 
effective  span.  This  can  be  done  by  Introducing  Intermediate  supports.  In 
this  study  such  supports  are  applied  to  the  slab  but  not  to  the  walls.  Such 
supports  are  referred  to  as  "expedient  upgrading”  and  may  consist  of  timber, 
steel ,  and  masonry. 

The  objective  here  Is  not  to  evaluate  the  particular  supports  used,  but 
rather  to  determine  the  reliability  of  the  shelter  when  an  Intermediate  support 
Is  provided.  Other  studies  have  been  devoted  to  the  design  and  experimental 
evaluation  of  expedient  upgrading  schemes  (Ref  1,2). 

The  first  floor  slab  was  designed  (Ref  3)  as  a  one-way  system  for  live 
loads  In  the  range  from  50  psf  to  250  psf.  This  represents  a  fairly  wide 
range  of  use  classes.  The  lower  bound  applies  to  classrooms  and  public  rooms. 


while  the  upper  bound  applies  to  Industrial  buildings,  e.g.,  light  manufacturing 
and  some  small  warehouses.  A  total  of  twelve  separate  cases  representing 
three  different  span  lengths  (12,  16,  and  20  ft)  and  four  different  design 
live  loads  (50,  80,  125,  and  250  psf)  were  considered.  Each  of  the  twelve 
basements  was  analyzed  as  expediently  upgraded  using  four  different  upgrading 
schemes.  This  resulted  In  sixty  shelters  which  Include  the  conventional  (un¬ 
upgraded)  design  as  the  base  case.  As  used  In  this  study,  an  expedient  up¬ 
grading  scheme  Involves  supporting  the  first  floor  slab  and  blocking  off  all 
openings  Into  the  basement. 

The  second  category  of  shelter  considered  Includes  four  conventional 
wood  frame  residences  with  basements.  These  are  real  buildings  whose  plans 
were  obtained  from  local  engineer/architect  offices.  Expedient  upgrading 
schemes  considered  In  this  portion  of  the  study  Include  the  "studwall"  and 
"post  and  beam."  The  objective  Is  to  reduce  the  effective  span  length  of  the 
joist  floor  over  the  basement.  Six  shelters  were  evaluated.  First,  each 
basement  was  evaluated  as  upgraded  using  the  studwall  scheme.  Second,  two  of 
the  basements  were  reevaluated  using  the  post  and  beam  concept.  The  process 
was  repeated  by  assuming  that  1  ft  of  soil  would  be  placed  over  the  first 
floor  for  radiation  protection.  Placing  2  ft  of  soil  would  significantly 
affect  the  strength  of  the  floor  system.  The  case  Involving  2  ft  of  soil  was, 
therefore,  not  considered. 

A  probability  of  survival  function  was  developed  for  each  shelter  and  each 
particular  upgrading  scheme.  The  method  used  In  determining  the  probability 
of  people  survival  Is  described. 

The  analysis  procedure  fbrmulated  and  used  In  this  study  consists  of  two 
parts.  The  first  part  Is  a  probabilistic  structural  analysis  which  determines 
the  probability  of  shelter  failure  (collapse).  This  analysis  Is  capable  of 
considering  all  structural  components  and  the  respective  failure  modes  of  each 
component.  For  example.  In  the  case  of  the  reinforced  concrete  slab,  both 
flexure  and  shear  are  considered  as  contributing  to  collapse.  The  probability 
of  failure  for  each  mode  acting  Independent  of  the  others  Is  determined  first. 
Correlation  between  them  Is  not  evaluated.  The  results  are  then  used  to 
determine  the  upper  and  lower  bounds  on  the  probability  of  failure  fbr  each 
component  and  then  for  the  structure  as  a  whole.  As  an  example,  see  the  analysis 
presented  In  Appendix  A. 
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The  second  part  Is  a  probabilistic  people  survival  analysis  which  Makes 
use  of  the  probability  of  structural  failure  results.  Casualty  mechanises 
considered  Include  debris  from  the  collapse  of  the  shelter  and  primary  blast. 
Probability  of  survival  against  primary  blast  Is  determined  on  the  basis  of 
available  casualty  data  (Ref  4).  This  report  Is  arranged  as  follows. 

Chapter  2  Includes  a  detailed  description  of  the  structural  analysis 
used  In  predicting  the  respose  of  reinforced  concrete  slabs  when  subjected  to 
blast  loading.  The  corresponding  probabilistic  analysis  Is  presented  In 
Chapter  3.  These  two  chapters  form  the  basis  of  a  computer  program  which  was 
started  In  the  previous  effort  (Ref  5)  and  then  modified  and  extended  In  the 
course  of  the  study  reported.  This  computer  program  Is  capable  of  computing 
the  probability  of  survival  of  people  located  In  basement  shelters  when  sub¬ 
jected  to  blast  produced  by  the  detonation  of  a  nuclear  weapon.  The  procedure 
used  Is  as  described  In  the  previous  paragraphs  of  this  chapter. 

In  Its  present  form  this  computer  program  can  analyze  basement  shelters 
In  which  the  roof  slab  Is  the  primary  structural  component  and  the  walls  are 
not  considered  In  the  analysis.  It,  therefore,  applies  to  cases  In  which  the 
walls  are  not  exposed  to  the  blast,  or  by  virtue  of  their  design  and  location 
are  substantially  stronger  than  the  slab.  The  program  consists  of  two 
separate  parts,  which  treat  the  following  problems: 

(1)  Basements  with  two-way  roof  slabs  and  with  membrane  resistance 

along  two  or  four  opposite  edges.  In  addition  to  membrane 

resistance,  four  support  conditions  may  be  considered. 

(a)  All  edges  simply  supported 

(b)  All  edges  fixed  (clamped) 

(c)  Long  edges  simply  supported,  short  edges  fixed 

(d)  Short  edges  simply  supported,  long  edges  fixed. 

(2)  Basements  with  two-way  roof  slabs,  without  membrane  resistance. 

These  four  support  conditions  may  be  considered. 

(a)  All  edge*  simply  supported 

(b)  All  edges  fixed 

(c)  Long  edges  simply  supported,  short  edges  fixed 

(d)  Short  edges  simply  supported,  long  edges  fixed. 
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The  computer  program  computes  the  probability  of  survival  for  people 
located  In  basements  when  subjected  to  blast  effects  produced  by  the  detonation 
of  a  nuclear  weapon  In  Its  Mach  region.  Megaton  or  klloton  weapons  may  be 
specified.  The  computer  program  has  some  of  the  following  features. 

(1)  It  predicts  the  upper  and  lower  bounds  on  the  probability 
of  component  collapse.  In  doing  this  both  the  flexural 
and  shear  modes  of  failure  are  considered. 

(2)  Predicts  the  probability  of  people  survival  based  on: 

(a)  Debris  effects  from  the  collapse  of  the  overhead  slab 

(b)  Blast  pressures  due  to  primary  blast. 

Slab  collapse  modes  on  which  the  debris  effects  are  based 
were  estimated  based  on  review  of  experimental  data. 

(3)  Considers  statistical  variability  In  the  following  parameters: 

Blast  load  parameters  -  ,  trf 

Structure  parameters  -  Ae,  A',  f',  f  ,  d,  d',  <t> 

5  5  c  y 

where  F1  »  peak  overpressure 

td  *  positive  phase  duration 

A„  ■  tension  steel 
s 

Aj  *  compression  steel 

f*  *  compressive  concrete  strength 

fy  ■  reinforcement  yield  strength 

4  *  undercapacity  factor  In  bending. 

It  Is  our  considered  opinion  that  this  computer  program  Is  superior  to 
any  that  exist  In  related  areas.  The  reasons  are: 

(1)  The  program  analyzes  actual  structures  and  makes  predictions 
on  the  basis  of  analytic  results.  No  scaling  Is  Involved. 

(2)  Parameter  variability  Is  considered  In  more  detail  and  on 
a  larger  scale  then  other  methods  (such  as  the  FAST  code. 

Reference  19,  for  example).  The  results,  therefore,  are 
more  reliable. 

(3)  The  program  Is  capable  of  evaluating  the  effectiveness  of 
different  expedient  upgrading  schemes  on  people  surviv¬ 
ability. 
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(4)  The  program  does  not  use  a  simulation  approach,  such  as 
Monte  Carlo,  for  example,  and  Is,  therefore,  quick  and 
economical  In  computer  usage. 

(5)  The  program  Is  oriented  specifically  to  the  civil  defense 
(national  security)  problem. 

This  computer  program  can  and  should  be  extended  to  Include  a  more  com¬ 
plete  set  of  structural  components,  l.e.,  walls,  columns,  and  girders.  This 
would  extend  Its  applicability  to  a  wider  class  of  structures  and  would  thus 
Increase  Its  utility. 

For  the  sake  of  clarity  and  generality  In  presentation,  the  structural 
analysis  procedure  given  In  Chapter  2  and  the  probabilistic  analysis  pro¬ 
cedure  given  In  Chapter  3  are  explained  with  reference  to  a  square,  two-way 
slab  fixed  along  the  edges. 

Reinforced  concrete  shelters  considered  here  are  described  In  Chapter  4 
which  also  Includes  a  description  of  the  expedient  upgrading  options  con¬ 
sidered.  People  survivability  results  are  summarized  In  Chapter  5.  Conclu¬ 
sions  and  recommendations  are  presented  In  Chapter  5  together  with  a  short 
summary  of  this  study. 

Analysis  of  residential  basements  Is  presented  In  Appendices  A  and  B. 
Appendix  C  contains  detailed  probability  of  failure  and  probability  of  sur¬ 
vival  results  for  the  reinforced  concrete  basement  shelters. 
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2.  STRUCTURAL  ANALYSIS 


2.1  STRUCTURE  IDEALIZATION  AND  LOADING 

The  general  form  of  the  structural  analysis  performed  In  this  study  Is 
described  In  Reference  6  and  discussed  In  relation  to  the  structure  shown 
In  Figure  1.  This  Is  a  portion  of  a  conventional  basement  which  serves  as 
a  personnel  shelter  against  the  effects  of  blast  produced  by  a  nuclear 
weapon  detonated  near  the  ground  surface.  With  the  entranceways  and  other 
openings  Into  the  basement  blocked  off.  the  structural  component  of  primary 
Interest  Is  the  first  floor  slab.  Its  collapse  will  result  In  casualties 
due  to  slab  debris  Impact  and  due  to  blast  pressures  and  velocities  pene¬ 
trating  basement  areas  where  people  would  be  located.  Basement  walls,  both 
Interior  and  exterior  (peripheral),  are  not  expected  to  fall  prior  to  the 
failure  of  the  first  floor  slab  and  are,  therefore,  not  considered  In  the 
analysis. 

The  roof  (first  floor)  slab  Is  modeled  as  a  single  degree  of  freedom 
system  whose  resistance  Is  a  piecewise  linear  function.  The  point  at  which 
response  Is  sought  Is  at  the  center  of  the  slab.  We  are  Interested  In  Its 
peak  deflection  when  the  slab  Is  subjected  to  a  time  dependent  load  over  Its 
surface.  We  are  also  Interested  In  the  peak  dynamic  reactions  distributed 
along  the  edges  of  the  slab.  The  blast  load  Is  approximated  using  the 
following  expression  (Ref  7): 


F(f) 


F1  N 


(1) 


where  F^  ■  peak  load  magnitude 

t^  ■  positive  phase  duration  of  the  blast  load. 

The  spatial  distribution  of  the  blast  load  Is  assumed  to  be  uniform  over  the 
surface  of  the  slab.  The  Interaction  of  the  blast  wave,  with  the  building 
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Basement  shelter*  plan  and  elevation  views 


above  the  basement  Is  assumed  not  to  modify  the  free- field  character  of  the 
blast  wave  to  any  significant  extent.  Therefore,  equation  (1)  need  not  be 
modified  for  this  effect. 

2.2  RESISTANCE  FUNCTION 

A  resistance  function  for  uniformly  loaded,  two-way  reinforced  concrete 
slabs  fixed  along  the  edges  Is  shown  In  Figure  2.  The  maximum  resistance 
In  the  elastic  range,  Rj  Is  assumed  to  be  developed  when  the  most  highly 
stressed  section  reaches  Its  plastic  resistance.  For  slabs  fixed  along  the 
edges  this  section  Is  along  the  long  edges.  For  a  square,  clamped  slab  for 
example,  R-j  Is  (Ref  6) 


R-|  -  29.2  M*b  (2) 


where  H®h  ■  negative  ultimate  moment  capacity  per  unit  width 
UD  at  the  center  of  the  long  edge. 

The  maximum  resistance  In  the  elasto-plastlc  range,  R2,  Is  determined  on 
the  assumption  that  the  ultimate  bending  moment  Is  developed  along  all  yield 
lines  representing  a  minimum  load  yield  pattern.  Thus  for  a  square  clamped 
slab  (Ref  6) 

*2  ■  T  (Mul  +  *u2  ♦  "u3  ♦  1u4>  <’> 


where  M  .  ■  total  ultimate  moment  capacity  along  mlspan 
Ul  section  parallel  to  edge  "a" 

M  -  *  total  negative  ultimate  moment  capacity 
along  edge  "a" 

M  ,  »  total  ultimate  moment  capacity  along  midspan 
UJ  section  parallel  to  edge  "b" 

M  ,  ■  total  negative  ultimate  moment  capacity  along  edge  NbM 
'**  (see  Figure  3  for  distribution  of  ultimate  moments) 

a  -  span  length  In  the  short  direction. 
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The  Midpoint  deflection  y1  Is 


"i 


where  k. ,  the  slab  stiffness  In  the  elastic  range  for  a  square,  clamped 
slab  for  example,  Is  (Ref  6) 


k 


1 


810  E  I, 

A 


(5) 


where  E  ■  modulus  of  elasticity  of  concrete 

I  -  moment  of  Inertia  of  a  unit  width  of  slab. 
0 


E  -  33w3/Z  /7T 
c 


(6) 


where  w  ■  unit  weight  of  concrete,  l.b/cu  ft 

f*  ■  ultimate  compressive  strength  of  concrete,  psl. 

I#  -  (5.5  pa  +  0.083) 


(7) 


where  p  *  average  reinforcement  ratio,  for  a  slab  with 
*  uniform  reinforcement,  p.  •  As/bd, 

where  As  Is  the  steel  area  per  unit  width,  b, 
and  d  Is  the  effective  depth  of  the  section. 


y2  *  yl 


(8) 


where  k„  the  slab  stiffness  In  the  elasto-plastlc  range  also  for  a  square, 
2  clamped  slab.  Is  (Ref  6) 


252  E  I 

4,-jr 
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Based  on  an  examination  of  experimental  results.  Reference  8  recommends 
that  the  failure  (Incipient  collapse)  deflection,  ym»  be  computed  as 

ym  ■  0.15a  (10) 


2.3  DYNAMIC  REACTIONS 

The  dynamic  reactions  along  the  edges  of  two-way  slabs  are  determined 
on  the  basis  of  the  assumption  that  the  distribution  of  the  Inertial  forces 
Is  the  same  as  the  assumed  deflected  shape  of  the  slab  and  the  resistance 
Is  uniformly  distributed. 

For  a  fixed,  uniformly  loaded  two-way  slab  with  a/b  =0.5  for  example, 
the  total  dynamic  reaction,  Va,  at  edge  "a"  in  the  elastic  range  is  (Ref  6) 

a 

V  «  0.05P  +  0.08R  (11) 

9 


The  corresponding  dynamic  reaction  at  edge  "b"  Is 

=  0.12P  +  0.25R 

D 

In  the  elasto-plastlc  range 

V,  «  0.04P  +  0.09R 

9 

Vb  =  0.09P  +  0.28R 

In  the  plastic  range 

=  0.04P  +  0.08R2 
Vb  =  0.11P  +  0.27R2 

where  P  *  abF^(t) 

R  *  R(y) 


(12) 


(13) 

(14) 


(15) 

(16) 

(17) 

(18) 
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It  was  assumed  In  this  study  that  the  dynamic  reactions  are  uniformly 
distributed  along  the  respective  edges  (Figure  4).  The  critical  shear  stress, 
7y,  was  computed  at  a  section,  d/2,  from  the  face  of  support  using  the 
approximation  shown  In  Figure  5.  Thus,  the  critical  shear  stress  along  edge 
"b"  Is 


where 


ua 


1 

** 

1 

<0 

d) 

(19) 

vua  (b  -  t.  - 

d) 

(20) 

~T~  w 

vub  *  H 

(21) 

V. 

V  a  — T- 

ua  aa 

(22) 

t  ■  support  (wall)  thickness 
d  »  effective  depth  of  slab. 


2.4  SHEAR  RESISTANCE 

The  shear  resistance  provided  by  the  concrete  can  be  computed  using  the 
following  expression  (Refs  9,10) 


'  (’• 


2500p„ 


Vud 


) 


bd 


(23) 


but  not  greater  than  3.5  bd  /T£t 
where  pw  ■  Aj/bd 

b  ■  width  of  section 

Vu  *  the  shear  at  the  section 

M  *  the  bending  moment  at  the  section  occurring 
u  simultaneously  with  Vu. 

The  quantity  w\  Is  not  to  he  taken  greater  than  1.0  In  computing  Vc< 
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The  shear  resistance  can  also  be  approximated  (Ref  10)  by 


Vc-2/TJbd  (24) 

According  to  Reference  9,  equation  (23)  Is  too  conservative  for  predicting 
structural  failure  due  to  shear  and  recommends  that  values  obtained  from  equa¬ 
tion  (23)  be  Increased  by  50%  when  used  for  this  purpose.  Consequently,  the 
following  expression  was  used  In  Investigating  shear  failure  of  slabs: 

'.■'•'h  <25> 

where  vm  ■  the  unit  capacity  of  the  slab. 

In  computing  vm,  the  ultimate  dynamic  strength  of  concrete,  f'  ,  was  used. 


1 .25  f ' 

C 


(26) 


2.5  DYNAMIC  RESPONSE  ANALYSIS 

Since  both  the  loading  and  resistance  are  complex  functions.  It  was 
necessary  to  use  a  numerical  procedure  to  obtain  the  peak  midpoint  deflec¬ 
tion  of  the  slab  and  the  peak  dynamic  reactions.  The  equation  solved  Is 

Kl„  Mt?  +  R(y)  -  F(t)  (27) 


where  *  the  load-mass  factor  which  Is  used  to  transform  the 
real  system  to  an  equivalent  single  degree  of  freedom 
system.  For  a  square,  clamped  reinforced  concrete 
slab,  Ky^  has  the  following  values  (Ref  6): 

Elastic  range,  Ky^  *>  0.63 
Elasto-plastlc  range,  Ky|  ■  0.67 
Plastic  range,  Kjjj  ■  0.51 
In  the  membrane  range,  KLM  was  taken  as  1.0 

Mj.  *  the  total  mass  of  the  slab 

R(y)  »  resistance 

F(t)  ■  load-time  history,  see  equation  (1). 


3.  PROBABILISTIC  ANALYSIS 

The  survivability  or  vulnerability  of  a  structure  to  a  given  load  Is  a 
matter  of  available  resistance  relative  to  the  Imposed  load.  If  the  load 
and  resistance  could  be  specified  exactly,  there  would  be  no  question  about 
predicting  survivability.  However,  due  to  uncertainties,  neither  the  load 
nor  the  resistance  can  be  specified  precisely  and  for  this  reason  surviv¬ 
ability  needs  to  be  expressed  In  terms  of  a  probability. 

For  a  structure  with  resistance  R  and  load  S,  where  R  and  S  are  random 
variables,  survival  Is  the  event  R>S  and  conversely,  failure  Is  the  event  R<Si 
If  fs  and  fR  are  respectively  the  probability  density  functions  of  applied 
load  and  resistance,  then  the  probability  of  failure,  P(F),  may  be  related  to 
the  overlapping  region  between  f  and  fD  (see  Figure  6).  Accordingly,  the 
probability  of  failure  Is  a  function  of  the  relative  position  between  vs  and 
u  (see  Figure  6)  where  y  and  y  are  the  expected  values  of  S  and  R,  res- 

I  5  i 

pectlvely.  The  probability  of  failure  also  depends  on  the  degree  of  uncer¬ 
tainty  (dispersion)  In  R  and  S  as  shown  In  Figure  7. 

Uncertainties  arise  due  to  variability  In  each  of  the  load  and  resistance 
parameters,  and  due  to  Imperfections  In  the  analytic  models  used  In  calculating 
load  and  resistance. 

3.1  PROBABILISTIC  ANALYSIS  APPROACH 

The  previous  discussion  points  out  the  importance  of  treating  the  problem 
of  survivability  evaluation  In  probabilistic  terms.  The  corresponding  general 
framework  for  doing  this  Is  described. 

For  a  given  structure.  Its  performance  function  Z  can  be  defined  as 

Z  -  R  -  S  (28) 


15 


and  f  f  and  f 


R  and  S  are  respectively  functions  of  several  variables  and,  therefore,  also 
Z,  thus 


Z  =  g  (X1t  X2 . XN)  (29) 

The  performance  limit,  l.e.,  the  minimum  level  of  performance  that  Is  required 

for  survival,  can  be  set  at  Z  *  z  .  If  z^  »  0,  then  Z  5  z„  defines  a  failure 

oo  o 

state  and  Z  >  zq  defines  a  survival  state.  Accordingly,  the  probability  of 
survival,  P(7)  Is 


P(7)  *  P(Z  >  zQ) 


(30) 


and  the  probability  of  failure,  P(F),  Is 

P(F)  =  1  -  P( F)  »  P(Z  <  zq)  (31) 


If  the  probability  distribution  (probability  density  function)  of  Z  is 
known,  then  the  probability  of  failure  would  be 


P(F) 


fz(z)dz  =  Fz(zQ) 


(32) 


Consider  the  case  In  which  R  and  S  are  known  to  be  lognormal  with  mean 
(expected)  values  yR  and  ys  and  coefficients  of  variation  nR  and  Assume 
also  that 

(1)  R  and  S  are  statistically  Independent 

(2)  The  probability  distributions  on  R  and  S  are  known 

(3)  The  uncertainty  Is  due  only  to  randomness  as  described  In 
the  distribution  fz(z) 

(4)  The  performance  limit  zQ  ■  0. 
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For  this  case,  the  performance  function  Is 


Z  -  *n(R/S)  (33) 

Since  R  and  S  are  lognormal  variates,  their  logarithms  are  normal  variates  and 
thus  Z  *  £n( R/S)  =  ZnR  -  ZnS  Is  also  a  normal  variate.  Therefore,  the  prob¬ 
ability  of  failure  Is 

P(F)-Fz(2o).t(^_L (34) 

where  *  (  )  *  cumulative  density  function  of  the  normal  distribution 
u2  *  mean  (expected)  value  of  Z 
02  ■  standard  deviation  of  Z. 

The  expected  value  and  standard  deviation  of  Z  are  obtained  as  (Ref  11) 

uz  »  jtnug  -  TjAn  (1  +  QR)J  -  jtnuj  -  ^£n  (1  + 

where  yR,  *  mean  values  of  R  and  S,  respectively 

flp,  ■  coefficients  of  variation  of  R  and  S,  respectively. 

Rearranging  R  and  S,  uz  becomes 


yZ 


■  Zn 


“s  V’T‘? 


a\  -  V(Z)  -  V(ZnR)  +  V(ZnS)  ■  Zn(1  +  fljj)  +  Zn(l  +  fl|) 
where  V(  )  *  variance  of  respective  parameter 

■ - r - n - a-T 


(36) 


(37) 


Substituting  pz  and  az  Into  equation  (34),  P(F)  becomes 


P(F)  *  * 


-Jtn 


^rV1  +4 

hjrr%\ 


fiT'o  +  «r)o  +  n^)1 


If  CJjj  and  are  not  too  large,  less  than  about  0.30,  then 

■tn  (vR/vs) 


P(F)  *  * 


(39) 


(40) 


The  task  of  determining  fz(z)  and  thus  Fz(z)  Is  generally  very  Involved. 
The  required  distribution,  fz(z).  Including  Its  parameters,  needs  to  be  de¬ 


rived  from  X| ,  Xz»  ..., 


XR  consistent  with  equation  (29).  Given  the  density 


functions  fx  (xj),  fx  (x?),  ....  fx  (xn),  the  cumulative  density  function 
(CDF)  of  Z  for  Independent  X-j ,  Xg, 


...»  Xp  would  be  - 


F,(z)  -  /  f v  f  ...  f  dx,  dx?  ...  dxn 

Z  X1  x2  xn  1  c 

{9(^1 1  X£ ,  ...,  Xjj  -  z} 


(41) 


from  which  the  probability  of  failure  would  be  obtained. 

The  effort  to  derive  the  exact  distribution  of  Z  through  equation  (41) 

Is  clearly  laborious  and  In  most  cases  Impractical,  because  the  density 

functions  f  ,  f  ,  ....  f  are  usually  not  known.  The  Information  on  these 
X1  x2  xn 

variables  Is  generally  limited  to  the  mean  values  and  coefficients  of  varia¬ 
tion.  However,  the  necessary  type  of  distribution  for  Fz(z)  may  be  prescribed, 
taking  Into  account  relevant  physical  factors  that  could  contribute  to  the 
distribution  form  with  consideration  for  mathematical  tractablllty. 
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With  regard  to  the  required  distribution  for  Z,  we  observe  the  following 
(Ref  12): 

(1)  When  P(F)  >  10“^,  the  calculated  P(F)  Is  approximately 
the  same  Irrespective  of  the  assumed  distribution  for  Z. 

(2)  When  P(F)  <  10"®,  the  calculated  P(F)  could  be  very 
sensitive  to  the  assumed  distribution  form. 

In  the  light  of  these  observations,  the  correct  probability  of  failure 

-3 

may  be  estimated  using  any  reasonable  distribution  for  Z  when  P(F)  >10  ; 

-5 

whereas  for  cases  where  P(F)  <  10  ,  a  correct  distribution  for  Z  would  be 
necessary  to  estimate  the  true  risk. 

In  the  case  where  R  and  S  are  not  Independent,  equation  (30)  would 
Involve  the  covariance  between  R  and  S  and  the  method  given  In  Reference  13 
would  need  to  be  used. 

3.2  ESTIMATING  MEANS,  VARIANCES,  AND  COVARIANCES 

For  the  general  function  of  a  random  variable  x,  l.e.,  y  ■  f(x),  the 
mean  and  variance  are 


E[f(x)] 


Jf f(x)g( 


x)dx 


(42) 


V[f(x)] 


Cf(x)  -  v]2g(x)dx 


(43) 


where  y  ■  E[f(x)] 

g(x)  ■  the  probability  density  function. 

As  discussed  earlier  In  this  narrative.  In  many  practical  applications 
g(x)  may  not  be  known  and  Information  may  be  limited  to  the  mean  and  variance 
of  the  original  variate  x.  Even  If  g(x)  Is  known,  the  Integrations  Indicated 
by  equations  (42)  and  (43)  may  be  difficult  to  perform.  For  these  reasons, 
approximate  expressions  fbr  the  mean,  variance,  and  covariance  may  be  obtained 
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by  expanding  known  functions  In  Taylor  series  and  neglecting  all  terms  except 
the  linear  terms.  If  R  Is  a  function  of  n  variables,  l.e.. 


R  *  f(x^ ,  *2*  •••» 


then  It  can  be  shown  (Ref  14)  that  the  expected  value  of  R  Is  - 


E(R)  *  f(x^ i  *2*  •••*  x^) 


where  x^  »  the  mean  value. 

Also,  the  variance  of  R,  V(R)  Is 


.  f  (aS_Y  V(x  )  ♦  2  V  f  >131 

M9V  '  i-l  k  3X1  8xj 


C  (x. ,  x.) 


where  C(x^,  x^)  *  covariance  between  x^  and  Xj. 

If  two  functions,  R-j  and  Rg,  are  functions  of  the  same  n  variables,  l.e.. 


Ri  ’  Vxr  x2-  •••’  xn) 


R2  ’  f2(xl*  x2 . xn> 


then  the  covariance  between  Rj  and  Rg  has  the  following  form  (Ref  14): 


n  3R,  3R, 

C(R,  •  R2>  ■  £  1J7  j4V(xl' 


3.3  MULTIPLE  FAILURE  MOOES 


If  there  Is  more  than  one  mode  of  ft 1 lure  for  a  given  structural 
component,  the  probability  of  survival  will  be  a  function  of  the  respective 
ftllure  modes.  Denoting  by  Rj  the  resistance  In  mode  1  and  by  S.|  the 
corresponding  load  effect,  the  survival  probability  Is  theoretically  (Ref  15) 

P(F)  -  P(R1  >  S1  n  R2  >  S2  n...n  Rn  >  Sn)  (49) 

Equation  (49)  applies  to  cases  where  the  failure  modes  are  Independent. 

When  the  failure  modes  are  highly  correlated,  the  probability  of  survival 
becomes 


P(7)  -  min  P(Rt  >  Sj)  (50) 

which  means  that  the  survival  of  the  component  Is  determined  by  the  weakest 
mode. 

3.4  ANALYSIS  OF  SYSTEM  SURVIVABILITY 

The  probability  of  survival  described  thus  far  refers  to  that  of  a 
single  structural  component.  The  survivability  of  a  complete  structure 
consisting  of  combinations  of  such  components  will  depend  on  the  surviv¬ 
abilities  of  these  components  and  the  manner  In  which  they  are  arranged 
and  connected.  The  analysis  will  need  to  consider  the  correlation  between 
the  components  and  the  degree  of  redundancy  of  the  system.  A  procedure 
for  accomplishing  this  Is  described  In  Reference  16. 

3.5  MEAN  ANO  VARIANCE  OF  THE  PEAK  MIDPOINT  DEFLECTION  OF  THE  SLAB 

This  section  contains  expressions  of  the  expected  value,  yp,  and  the 
variance,  V(yp)  of  the  peak  midpoint  deflection,  yp,  In  the  respective 
ranges  of  response,  l.e.,  elastic,  elasto-plastlc,  and  plastic. 
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When  the  magnitude  of  the  blast  load  Is  such  that  the  response  of  the 
slab  Is  In  the  elastic  range,  then  (See  Fig.  2) 

R(y)  *  yRi/Yi  <51) 

In  accordance  with  equations  (27)  and  (45)  the  expected  value  of  the  peak 
midpoint  deflection  can  be  expressed  as  a  function  of  the  following  parameters, 
i  • 


yp  *  ff^j »  td,  R1t  y-j ,  Mt)  (52) 

As  Indicated  previously,  yp  was  determined  numerically.  The  variance  of  yp 
was  also  determined  numerically  using  the  following  expression: 


C(RrMt}  +  3^  c(yi  ’V 


Similarly,  when  the  response  Is  In  the  elasto-plastlc  range, 

R(y)  *  Rj  +  y*-_— ~  (y  “  y {54) 
Yp  *  f (^i »t(j,Rl ,^2'Yl ^2’^)  (55) 
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3y_  3y  3y  3y  3y_  3y 

+  3^  W;  C(R1  »Mt)  +  3^  3yf  C(R2*y1 }  +  3Rj  3yJ  C(R2’y2} 

3y  3y  3y  3y  3y  3y 

+  ^  3H^  C(R2.Ht)  +  gj*  Cfy,  ,y2)  +  ^  C[y,  ,Mt) 


When  the  response  Is  In  the  plastic  range,  R(y)  ■  R?  and 


>p  *  ^VW^'VVVV  <57> 


3y_  3y_  3y_  3y_  3y_  3y_ 

+  wfXC("1,"t>  *^^C<RZ,,1)  *^1^C<R2,y2) 

3y„  3y_  3y_  3y_  3jL  3y„ 

♦a^aiJ^Vt* ’ijfuJ  c{yi  •y2)  *  ?jf  aij  c(yi  •"»> 

9ym  i 

+  C(y2*"t)j 
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3.6  MEAN  AND  VARIANCE  OF  THE  MAXIMUM  SHEAR  STRESS 


This  section  contains  expressions  of  the  expected  value,  and  the 
variance,  V(vup)  of  the  maximum  shear  stress  In  the  respective  ranges  of 
slab  response,  l.e.,  elastic,  elasto-plastlc,  and  plastic. 

As  Indicated  previously  (see  Section  2.3),  the  shear  stress  In  the 
concrete  along  a  given  edge  of  a  slat.  Is  computed  using  equations  (19)  or 
(20).  When  the  magnitude  of  the  load  Is  such  that  response  Is  In  the  elastic 
range  then  the  peak  shear  stress  can  be  expressed  as  a  function  of  the 
following  parameters,  i.e., 

vUp  a  f^l*td*^l*yl,‘^ 

The  magnitude  of  was  determined  numerically  for  each  edge  and  the  maximum 
of  the  two  was  used  In  the  analysis.  The  variance  of  v  was  also  determined 

(j  | / 

numerically  using  the  following  expression: 


Similarly,  when  the  response  Is  In  the  elasto-plastlc  range, 

vUp  *  f(T^  ,t(j,l?i  ,T?2»y^  »y2»^ 
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When  the  response  Is  In  the  plastic  range. 


In  the  above  expressions  the  blast  load  parameters,  F-j  and  td,  are  Inde¬ 
pendent  of  the  resistance  function  parameters  R-|»  Rg.  y-j»  ym»  and  the 
mass  of  the  slab,  M^.  Although  t^  Is  a  function  of  ,  the  covariance  using 
equation  (48)  turned  out  to  be  zero.  The  ultimate  (failure)  deflection,  ym» 
see  equation  (10),  Is  Independent  of  the  other  resistance  function  parameters 
because  the  span  length  "aM  Is  taken  as  a  constant. 

In  the  above  expressions  for  the  variances  of  yp  and  vup,  the  constituent 
variances  and  covariances  were  determined  In  closed  form.  The  partial 
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derivatives  were  obtained  by  means  of  numerical  differentiation  using 


3yp .  (Vi+i-(Vt-i 

3x.j 


(65) 


where  *  (F,  .td.ym.Mt.Qt.Qft^ft) 

Ah  =  differentiation  Increment,  Ah  was  taken  as  0.10  of 
one  standard  deviation  for  each  of  the  variables. 


3.7  PROBABILITY  OF  SHELTER  FAILURE 

For  the  reinforced  concrete  basement  shelters  considered  In  this  study 
(see  Figure  1),  the  first  floor  (overhead)  slab  Is  considered  to  be  the 
primary  structural  component.  Its  collapse  will  result  In  casualties.  The 
probability  of  failure  of  the  slab  was  determined  on  the  basis  of  the  theory 
presented  In  Sections  3.1  to  3.3.  Specific  expressions  used  In  the  analysis 
are  given. 

3.7.1  Probability  of  Slab  Failure 

As  Indicated  In  Section  3.3,  In  the  case  of  two  failure  modes,  the 
probability  of  slab  failure,  P(F),  is 

P(F)  =  1  -  [1  -  P(Fb)][l  -  P(Fy)]  (66) 

If  the  modes  are  Independent,  and 

P(F)  *»  max  [P(Fb),  P(Fy)3  (67) 

If  the  modes  are  highly  correlated. 

In  equations  (66)  and  (67),  P(Fb)  ■  probability  of  failure  due  to 
bending  (flexure),  see  Section  3.7.2,  and  P(Fy)  *  probability  of  failure 
due  to  shear,  see  Section  3.7.2. 
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The  real  failure  probability  of  the  slab  Is  between  these  two  cases* 

1  *e  •  * 

max  [P(Fb),  P(Fy)]  <  P(F)  <  1  -  [1  -  P(Fb)][l  -  P(Fy)]  (68) 

In  the  analysis  performed  In  this  study,  both  bounds  were  calculated 
and  are  Included  with  the  results. 

3.7.2  Probabilities  of  Failure  Due  to  Bending  and  Shear 

In  accordance  with  equation  (39),  probabilities  of  failure  due  to 
flexure  and  shear  were  computed  using  the  following  expressions: 


where  y_  and  y  were  defined  previously,  see  Section  3.5 

m  p 

n  ■  coefficient  of  variation  of  y 
yP  p 

ft  *  coefficient  of  variation  of  y 
ym  m 


where  vffl  and  vup  were  defined  previously,  see  equation  (25)  and  Section  3.6 

flv  ■  coefficient  of  variation  of  v 
op 

ft  ■  coefficient  of  variation  of  vm. 


3.8  PROBABILITY  OF  PEOPLE  SURVIVAL 


Casualty  mechanisms  considered  In  this  study  Include  primary  air  blast 
and  debris  Impact  due  to  the  collapse  of  the  first  floor  slab. 

Body  damage  due  to  primary  air  blast  results  when  the  blast  wave 
engulfs  the  body.  In  such  a  case,  movement  of  different  tissue  masses  pro¬ 
duces  shear  waves  which  accelerate  different  organs  and  different  parts  of 
the  same  organ  to  different  velocities.  This  results  In  strains  and  fre¬ 
quently  In  ruptures.  Air-filled  organs  such  as  the  lungs  are  especially 
susceptible  to  this  type  of  damage.  Data  used  to  estimate  the  probability 
of  survival  from  this  effect  were  obtained  from  Reference  4  and  are  repro¬ 
duced  In  Table  1. 


TABLE  1.  PROBABILITY  OF  SURVIVAL  FROM  PRIMARY  BLAST,  P(SnK) 
(30  DAYS  AFTER  EXPOSURE)  pD 


Blast  Probability 

Overpressure  of  Survival  (<) 


40 

97.6 

50 

88.0 

60 

72.0 

70 

51.0 

80 

33.0 

100 

11.0 

120 

3.0 

150 

0.6 

Probability  of  survival  against  debris  due  to  the  collapse  of  the 
overhead  (first  floor)  slab  Is  determined  using  the  theorem  of  total  prob¬ 
abilities  as 


P(Ssc)  -  P(S|7)P(?)  +  P(S| F)P(F)  (71) 

where  P(S._)  »  probability  of  people  survival  against  structural  collapse 

P(S|T)  »  probability  of  people  survival  given  that  the  shelter 
does  not  fall  (collapse) 
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P(F)  *  probability  of  shelter  (structure)  survival 

P(S|F)  »  probability  of  people  survival  given  that  the  shelter 
falls  (collapses) 

P(F)  ■  probability  of  shelter  collapse  *  1  •  P(F). 

In  equation  (71)  the  probability  of  structural  failure,  P(F),  Is 
determined  as  described  previously.  P(S|T)  and  P(S|F)  are  determined  as 
described  In  the  following  paragraphs. 

Experimental  results  (Refs  17,18)  Indicate  that  two-way  reinforced  con¬ 
crete  slabs  falling  under  blast  loading  do  not  generally  break  up  Into 
separate  pieces.  The  concrete  cracks  throughout,  but  most  of  the  large 
pieces  remain  loosely  connected  to  the  reinforcing  steel.  On  the  other 
hand,  one-way  slabs  falling  under  blast  loading  do  not  remain  connected,  but 
rather  break  up  catastrophically  (Ref  2). 

At  overpressures  prior  to  the  total  collapse  of  a  slab  (one-way  or 
two-way),  spallation  Is  expected  to  occur  and  shelter  occupants  will  be  hit 
by  spalled  pieces  of  concrete.  Injuries  will,  therefore,  occur;  however, 
these  are  not  expected  to  be  at  the  fatality  level.  It  Is  assumed  that 
prior  to  slab  collapse  no  fatality  level  casualties  are  produced  by  slab 
debris  and  thus  P(S|T)  ■  1.0. 

When  people  are  uniformly  distributed  In  all  shelter  areas  and  are  lying 
down  at  the  time  of  the  attack,  then  P(S|F)  can  be  directly  related  to  the 
floor  area  unaffected  by  collapse.  In  other  words,  when  an  overhead  slab 
collapses  people  in  areas  affected  by  the  collapsed  debris  become  fatality 
level  casualties  while  people  In  areas  unaffected  by  collapsed  debris  are 
survivors  and/or  Injured. 

For  personnel  shelters  with  two-way  overhead  slabs,  a  procedure  for 
computing  P(S|F)  Is  given  In  Reference  5.  A  procedure  for  computing  P(S|F) 
for  shelters  with  one-way  overhead  slabs  Is  described  next. 

Figure  8  shows  an  expediently  upgraded  basement  shelter.  The  particular 
upgrading  consists  of  a  series  of  6  In.  x  10  In.  beams  and  columns  which  are 
used  to  support  the  existing  slab  such  that  four  smaller  slabs  are  produced. 
For  this  concept,  P(S|F)  Is  computed  as  follows: 
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(72) 


P(S|F)  -  Sf 


where  »  total,  clear  floor  area 

A  ■  floor  area  that  would  be  affected  by  debris  from 
a  the  collapse  of  the  slab. 

For  the  case  where  a  -  12  ft  and  b  ■  24  ft,  the  spacing  of  the  beams 
Is  such  that  a'  *  4.08  ft  and  b'  *  10.08  ft,  then 

Af  -  (12  -  0.83 ) ( 24  -  0.83)  =  258.81  sq  ft 

A  -  4(4.08) (1 0.08)  -  164.51  sq  ft 

a 

Sf  *  0.36. 

3.9  COMBINED  PROBABILITY  OF  SURVIVAL 

It  Is  reasonable  to  assume  that  an  Individual  Injured  from  several 
weapon  effects  would  have  a  lesser  probability  of  survival  then  If  the 
Injury  was  due  to  one  of  the  effects.  Information  on  how  the  simultaneous 
action  of  several  effects  from  a  single  weapon  combine  to  result  In  a  prob¬ 
ability  of  survival  Is  not  known  at  this  time.  For  this  reason,  survival 
probabilities  from  different  effects  are  combined  In  this  report  as  In¬ 
dependent  phenomena.  Thus  the  probability  of  people  survival,  P(S),  against 
primary  blast  and  slab  collapse  is  computed  as 


4.  DESCRIPTION  OF  SHELTERS 


4.1  BASIC  STRUCTURE 

Reinforced  concrete  shelters  considered  In  this  study  were  Illustrated 
In  Figure  1.  This  Is  a  portion  of  an  engineered  basement  of  a  low-rise 
reinforced  concrete  building.  The  first  floor  slab  Is  a  one-way  reinforced 
concrete  slab  with  Its  top  surface  at  grade.  The  slab  Is  simply  supported 
along  Interior  reinforced  concrete  walls.  Twelve  designs  were  performed 
(see  Table  2)  for  several  combinations  of  design  live  load  and  span  length. 
The  design  live  load  ranges  from  50  to  250  psf  and  the  span  (short  direction) 
from  12  ft  to  20  ft. 

4.2  EXPEDIENT  UPGRADING  SCHEMES 

Four  types  of  expedient  upgrading  schemes  were  considered  with  each  of 
the  twelve  slabs  given  In  Table  2.  As  Illustrated  In  Figure  9,  scheme  A  Is 
the  basic,  conventional  slab  and  schemes  B  through  E  represent  expedient 
upgrading  schemes  In  the  order  of  Increasing  strength.  Upgrading  Is  accom¬ 
plished  by  reducing  the  basic  slab  to  a  series  of  smaller  slabs.  This  would 
be  done  by  providing  supports  along  the  dash  lines  shown  In  Figure  9. 

Supports  that  may  be  used  for  this  purpose  were  shown  In  Figure  8. 

4.3  ANALYSIS  DATA 

Table  3  contains  structural  data  used  In  the  analysis  of  this  set  of 
shelters.  The  various  cases  and  expedient  upgrading  schemes  are  Identified 
In  Table  2  and  Figure  9. 

As  Indicated  earlier,  the  basic  slab  was  designed  as  a  one-way  slab. 

By  making  use  of  the  temperature  reinforcement  which  Is  placed  orthogonal 
to  the  main  reinforcement,  each  slab  In  each  expedient  upgrading  scheme  was 
analyzed  as  a  two-way  slab.  In  Figure  9,  A$3  Is  the  main  reinforcement  and 
A$i  Is  the  temperature  reinforcement. 
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TABLE  2.  STRUCTURAL  PARAMETERS  FOR  ONE-WAY  SLAB*  (REF  3) 


Yield  strength  of  reinforcing  steel,  fy  ■  60,000  psi. 
*$ee  Figure  1  for  slab  configuration. 


As  shown  In  Table  3,  the  effective  depth  of  slab,  reinforcing  steel 
areas,  ultimate  strengths  of  concrete  and  steel,  and  the  capacity  reduction 
factor  are  treated  as  random  variables.  The  other  parameters,  l.e.,  span 
lengths  and  the  slab  thickness,  are  treated  as  constants.  Coefficients  of 
variation  were  determined  based  on  data  given  In  Chapter  7  of  Reference  5. 
Coefficients  of  variation  of  peak  overpressure  and  load  duration,  see 
equation  (1),  were  taken  as  1.0%. 

In  Table  3,  parameters  "a"  and  "b"  are  center  to  center  of  supports 
dimensions,  while  "ac"  and  "bc"  are  clear  span  dimensions  after  the  place¬ 
ment  of  expedient  upgrading  supports.  Values  of  "ac"  and  "bc"  were  used  in 
computing  by  means  of  equation  (72). 

In  computing  for  the  basic  slab  (scheme  A,  Figure  9),  A^  and  Aa» 
see  equation  (72),  were  assumed  not  to  be  the  same.  In  computing  A  ,  it 
was  assumed  that  the  portion  of  the  floor  bounded  by  the  walls  and  a  line 
6  In.  from  the  walls  would  be  essentially  free  from  debris  effects.  Thus 
for  this  one  case  only 


Af  -  (b  -  0.83)(a  -  0.83) 
A  *  (b  -  1.83) (a  -  1.83) 

3 


(74) 
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REINFORCED  CONCRETE  SLAB  ANALYSIS  DATA 
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Ultimate  dynamic,  compressive  strength  of  concrete,  f£c  =  3.9063  ksi.  Coefficient  of  variation  of  f£c;  «  17.58% 
Ultimate  dynamic  yield  strength  of  steel,  fjy  *  84.018  ksi.  Coefficient  of  variation  of  fj y;  Gy  *  9.22%. 
Unit  weight  of  concrete,  w  *  150  pcf. 

Capacity  reduction  factor,  $  =  1.0.  Coefficient  of  variation  of  4;  fy,  •  4.7%. 


Ultimate  dynamic,  compressive  strength  of  concrete,  f£c  =  3.9063  ksl.  Coefficient  of  variation  of  f$c;  0^  *  17.58% 
Ultimate  dynamic  yield  strength  of  steel,  fdy  *  84.01R  ksl.  Coefficient  of  variation  of  fdy;  ny  -  9.22%. 
Unit  weight  of  concrete,  w  -  150  pcf. 

Capacity  reduction  factor.  >  =  1.0.  Coefficient  of  variation  of 


5.  PEOPLE  SURVIVABILITY  RESULTS 


This  chapter  contains  a  summary  of  results  on  the  probability  of 
survival  of  people  when  located  in  expediently  upgraded,  reinforced  concrete 
basements  and  subjected  to  the  blast  effects  of  a  1-MT  surface  burst. 

Casualty  mechanisms  include  debris  from  the  collapse  of  the  overhead  slab 
and  primary  blast.  Results  were  generated  using  a  computer  program  which 
was  developed  during  the  previous  study  (Ref  5)  and  modified  and  extended 
in  the  course  of  this  study.  Specific  shelters  considered  are  described 
In  Chapter  4. 

This  chapter  explains  how  the  results  were  produced.  Case  1 A  Is  used 
for  Illustration.  The  results  are  summarized  In  Table  4.  A  complete  set 
of  detailed  results  Is  included  In  Appendix  C. 

Case  1A  (see  Table  3)  Is  the  basic  conventional,  one-way  slab  over  the 
basement  designed  to  resist  Its  own  weight  plus  a  live  load  of  50  psf.  Its 
plan  dimensions  are  12  ft  by  24  ft  and  the  total  thickness  is  5.0  In.  When 
subjected  to  a  uniformly  distributed  blast  loading  over  Its  surface,  the 
possible  failure  modes  are  flexure  and  shear.  Probabilities  of  failure 
based  on  flexure  and  shear  when  acting  Independent  of  each  other  are  shown 
In  Figure  10.  Note  that  shear  and  flexure  are  both  Important  In  producing 
failure  and,  therefore,  both  need  to  be  considered.  However,  since  we  do 
not  know  how  these  two  failure  modes  correlate  then  the  best  that  can  be  done 
Is  to  bound  the  actual  failure  probability  as  discussed  In  Sections  3.3  and 
3.7.1.  Thus,  using  equations  (66)  and  (67),  the  upper  and  lower  bounds  on 
the  failure  probability  fbr  this  slab  were  computed  and  are  shown  In 
Figure  11.  Corresponding  bounds  on  the  probability  of  people  survival  are 
shown  In  Figure  12.  They  were  determined  using  equation  (71). 

It  Is  useful  to  compare  the  effectiveness  of  the  various  expedient  up> 
grading  schemes  on  people  survivability.  This  Is  done  for  Case  1  In  Table  5 
which  contains  overpressure  ranges  at  the  90%  and  50%  probability  of 
survival  levels. 
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TABLE  4.  SUMMARY  OF  RESULTS,  REINFORCED  CONCRETE  BASEMENT  SHELTERS  (continued) 
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SUMMARY  OF  RESULTS,  REINFORCED  CONCRETE  BASEMENT  SHELTERS  (continued) 


TABLE  4.  SUMMARY  OF  RESULTS,  REINFORCED  CONCRETE  BASEMENT  SHELTERS  (concluded) 
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OVERPRESSURE 


TABLE  5.  OVERPRESSURE  RANGES  AT  THE  90% 
AND  50%  PROBABILITY  OF  SURVIVAL  (CASE  1) 


Probability 

Expedient  Upqradlnq 

Scheme 

Level 

'  1 

- B  — 

D  E 

90% 

0.39-2.64 

0.55-4.18  0.87-11.73 

0.86-10.42  1.58-18.20 

50% 

2.77-6.60 

4.57-12.90  12.47-26.10 

11.21-  19.96- 

For  schemes  D  and  E  the  upper  bonds  are  not  Included  at  the  50%  prob¬ 
ability  of  survival  level.  The  reason  Is  that  the  curves  are  very  flat  after 
a  certain  point,  and  do  not  intersect  the  particular  probability  value  In 
the  overpressure  range  of  Interest  (see  Figure  12,  for  example).  This  Is 
also  the  reason  why  numbers  are  missing  In  certain  columns  In  Table  4.  It 
Is  evident  from  these  results  that  expedient  upgrading  can  be  very  effective 
in  providing  protection.  The  50%  survival  probability  at  2.77  psi  for  the 
basic  slab  becomes  almost  20  psi  when  upgrading  scheme  E  is  employed. 

Similar  trends  will  be  noted  for  the  other  slabs. 

If  It  is  a  matter  of  choosing  between  two  basements  for  expedient  up¬ 
grading,  then  obviously  the  one  that  was  designed  for  a  higher  live  load 
should  be  chosen  assuming  that  both  are  In  good  physical  condition  and  the 
design  loads  are  known  fbr  each.  The  key  item  In  expedient  upgrading  Is 
the  correct  design  of  the  supporting  system  and  Its  correct  Implementation. 

It  is  recommended  that  experimental  studies  be  Initiated  whose  objective 
would  be  to  generate  experimental  data  on  the  response  of  reinforced  concrete 
slabs  subjected  to  dynamic  loadings.  He  need  experimental  data  In  the  res¬ 
ponse  range  approaching  failure.  Available  experimental  data  on  the  shear 
failure  of  slabs,  the  distribution  of  reactive  forces  along  the  supports,  and 
the  Interaction  of  flexure  and  shear  prior  to  and  at  the  point  of  failure. 

Is  especially  limited  at  this  time.  Reliable  data  would  aid  In  the  develop¬ 
ment  of  accurate  failure  theories  and  also  In  the  development  of  design  and 
Implementation  criteria  for  expedient  upgrading  schemes. 

Since  the  upper  and  lower  bounds  on  the  failure  probability  are  fairly 
far  apart  for  most  cases  studied  (see  Figure  12),  It  becomes  useful  to 


determine  the  correlation  that  exists  between  the  two  failure  modes  and 
then  to  determine  the  actual  failure  probability.  This  was  not  done  In 
this  study  because  the  methodology  for  doing  this  was  not  available.  For 
the  present.  It  Is  believed  that  the  lower  bound  should  be  used  as  a  con¬ 
servative  measure  of  the  probability  of  people  survival. 


6.  SUMMARY.  CONCLUSIONS.  AND  RECOMMENDATIONS 


6.1  SUMMARY 

The  study  reported  here  was  concerned  with  predicting  the  probability 
of  survival  of  people  located  In  expediently  upgraded,  conventional  basements 
when  subjected  to  the  blast  effects  produced  by  the  detonation  of  a  1-MT 
weapon  near  the  ground  surface.  Two  categories  of  basements  are  considered, 
i.e.,  basements  of  engineered  buildings  and  basements  of  single  family 
residences. 

The  first  category  Includes  basements  of  low-rise  engineered  buildings. 
The  basement  walls  are  unexposed  and  the  first  floor  slab  (the  slab  over 
the  basement)  Is  at  grade.  The  first  floor  slab  Is  the  primary  structural 
component  for  the  basement  as  far  as  protection  from  the  blast  Is  concerned. 
Its  collapse  will  result  In  casualties  due  to  debris  Impact  and  blast  winds 
and  pressures  entering  shelter  areas  when  the  shelter  envelope  Is  breached. 

The  slab  over  the  basement  was  designed  as  a  one-way  reinforced  con¬ 
crete  system  for  live  loads  of  50,  80,  125,  and  250  psf  and  span  lengths  of 
12,  16,  and  20  ft.  This  Includes  a  total  of  twelve  conventional  basements 
having  a  representative  range  of  use  classes.  Each  of  the  twelve  basements 
was  analyzed  as  expediently  upgraded  using  four  different  upgrading  schemes. 
Upgrading  Is  accomplished  by  providing  supports  that  reduce  the  effective 
span  of  the  slab  and  by  blocking  off  all  openings  Into  the  basement.  This 
resulted  In  sixty  shelters  of  different  strengths  which  Include  the  con¬ 
ventional  ,  unupgraded  slab  as  the  base  case. 

The  second  category  Includes  four  conventional  wood  frame  residences 
with  full  basements.  These  basements  were  also  analyzed  as  expediently  up¬ 
graded.  Again,  upgrading  consisted  of  providing  Intermediate  supports  for 
the  joist  floor  system,  blocking  off  all  openings  Into  the  basement,  and 
mounding  the  building  with  soil  up  to  the  first  floor  level,  about  2  to  3  ft 
from  grade.  Expedient  upgrading  Included  "studwall"  and  "post  and  beam" 
concepts.  Six  shelters  were  analyzed.  First,  each  basement  was  analyzed 


as  upgraded  using  the  studwall  scheme;  second,  two  of  the  basements  were 
reanalyzed  using  the  “post  and  beam"  concept. 

A  probability  of  survival  function  was  developed  for  each  shelter  and 
each  upgrading  scheme.  The  procedure  used  to  accomplish  this  consists  of 
two  parts.  The  first  Is  a  probabilistic  structural  analysis  which  determines 
the  probability  of  failure  to  the  shelter  envelope.  The  second  part  Is  a 
probabilistic  people  survival  analysis  which  considers  two  casualty  producing 
mechanisms,  l.e.,  debris  Impact  and  associated  effects  from  the  collapse  of 
the  primary  structural  system  and  primary  blast.  The  probability  of  struc¬ 
tural  failure  Is  made  use  of  In  computing  the  probability  of  survival 
against  debris  effects. 

Th'  analysis  method  briefly  described  above  was  formulated  In  the  course 
of  this  study  and  the  previous  effort  (Ref  5),  and  a  portion  of  It  was  com¬ 
puterized.  The  computer  program  Is  capable  of  analyzing  reinforced  concrete 
shelters  of  the  type  discussed  previously  In  this  chapter,  and  of  predicting 
the  probability  of  survival  for  shelter  occupants  against  the  effects  of 
blast. 


6.2  CONCLUSIONS 

On  the  basis  of  the  study  conducted  and  results  obtained,  the  following 
conclusions  are  made: 

(1)  A  method  for  computing  the  probability  of  failure  of 
structures  (Ref  15)  was  examined  and  applied  to  the  analysis 
of  personnel  shelters.  This  method  Is  capable  of  condslder- 
Ing  all  structural  components  making  up  the  structure  and 
the  respective  failure  modes  of  each  component.  This  method 
Is  the  most  rational  that  Is  available  at  this  time  and  the 
results  are  believed  to  be  the  most  reliable  of  those  pro¬ 
duced  In  this  subject  area  to  date. 

(2)  Expedient  upgrading  can  be  very  effective  In  Increasing  the 
live  saving  potential  of  conventional  reinforced  concrete 
basements.  Conventional  basements  with  one-wey  reinforced 
concrete  overhead  floor  systems  designed  for  a  live  load  of 
50  psf  can  be  expediently  upgraded  to  result  In  a  probability 
of  survival  of  50%  at  20  psl.  Slabs  originally  designed  for 
a  higher  live  load,  125  psf  for  example,  can  be  easily  up¬ 
graded  to  result  In  a  probability  of  survival  of  50%  beyond 
30  psl.  It  can,  therefore,  be  supposed  that  two-way  floor 
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systems,  which  are  generally  stronger  than  one-way  floor 
system,  have  the  capacity  of  providing  even  greater  pro¬ 
tection  when  effectively  upgraded. 

(3)  Conventional  wood  framed  basements  are  capable  of  being 
expediently  upgraded  to  provide  protection  beyond  5  psi 
at  a  probability  of  survival  of  50%  or  better. 

(4)  The  analysis  procedure  presented  here  Is  sufficiently 
general  and  can  be  readily  extended  to  include  the  Influence 
of  other  hazards  which  accompany  a  nuclear  weapon  attack, 
l.e.,  prompt  nuclear  radiation,  fires,  and  fallout  radia¬ 
tion. 

(5)  A  capability  should  be  formulated  and  Included  In  the 
analysis  to  study  the  effects  of  evasive  action  taken  by 
shelter  occupants  on  the  probability  of  survival. 


6.3  RECOMMENDATIONS 

The  previous  effort  (Ref  5 )  and  the  study  reported  here, have  been  very 
useful  In  formulating  the  people  survivability  problem  on  a  rational  prob¬ 
abilistic  basis.  The  approach  Is  very  promising  and  if  allowed  to  develop 
further  will  produce  a  reliable  computational  tool  for  the  rating  of  shelter 
spaces,  evaluating  alternative  shelter  systems,  and  for  performing  damage 
limiting  studies.  With  this  end  In  mind,  two  tasks  are  recommended. 

6.3.1  Experimental  Task 

There  is  a  need  for  experimental  studies  to  be  Initiated  to  generate 
data  for  a  better  understanding  of  how  reinforced  concrete  structural  com¬ 
ponents  respond  In  the  range  approaching  failure,  l.e.,  what  failure  modes 
are  introduced,  how  they  Interact  with  each  other  for  different  slabs,  what 
Is  the  Influence  of  boundary  conditions  and  loadings  on  the  modes  of  failure. 
Experimental  studies  should  be  conducted  to  generate  data  capable  of  Im¬ 
proving  the  current  formulations  of  the  following  failure  criteria: 

•  Failure  criteria  for  horizontally  oriented  reinforced 
concrete  slabs;  one-way  and  two-way  floor  systems 

•  Failure  criteria  for  vertically  oriented  slabs  In 
contact  with  the  soil;  basement  walls 

•  Failure  criteria  for  columns  and  beams. 


A  test  plan,  outlining  the  number  of  experiments  that  would  be  needed 
to  produce  the  necessary  data  cannot  be  produced  at  this  time.  The  first 
task  would  be  a  review  of  all  available  experimental  data  on  this  subject. 
With  this  task  completed  It  would  be  possible  to  outline  a  preliminary  test 
plan. 


6.3.2  Analytic  Task 

The  computer  program  developed  thus  far  must  be  further  developed  to 
Include  the  capability  to  analyze  the  following  structural  systems  and  to 
Include  related  aspects: 

(a)  Individual  Structural  Components 

In  addition  to  the  library  of  structural  components  Included  In  the 
computer  program  at  this  time,  the  following  should  be  Implemented: 

Flat  slabs 
Flat  plates 
One-way  slabs 

Beams  (steel,  reinforced  concrete) 

Columns  (steel,  reinforced  concrete) 

Composite  steel  and  concrete  systems 
Masonry  systems. 

(b)  Weapon  Effects  Hazards 

In  addition  to  blast  and  debris  effects  Included  In  the  computer 
program  at  this  time,  other  nuclear  weapon  effects  and  Indirect  hazards 
should  be  Included. 

Prompt  nuclear  radiation 
Fallout  radiation 
Ground  shock 
Fires. 

(c)  Casual ty  Data 

Available  data  for  estimating  the  level  of  casualty  experienced  by 
Individuals  against  the  various  hazards  should  be  reviewed  with  the  object 
of  making  the  current  casualty  predicting  process  more  reliable. 

All  of  the  aspects  outlined  would  be  considered  within  the  probabilistic 
framework  outlined  In  Section  3, 
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APPENDIX  A 


PROBABILITY  OF  PEOPLE  SURVIVAL  IN  THE  BASEMENT 
OF  A  WOOD  FRAME  RESIDENCE 

Procedure  used  in  determining  the  probability  of  people  survival  In 
basements  of  single-family,  wood  frame  residences  Is  presented  In  this 
appendix.  The  particular  building  analyzed  Is  the  "Dunes"  house  (Ref  20). 

This  Is  a  one-story,  single-family  frame  residence  with  a  full  basement. 

The  floor  system  over  the  basement  Is  approximately  1  ft  above  grade. 

The  analysis  described  Is  for  a  basement  expediently  upgraded  by  pro- 
Idlng  stud  walls  In  the  basement  as  additional  supports  for  the  floor 
system,  blocking  windows  and  doors,  and  mounding  the  structure  on  the  outside 
up  to  the  first  floor  level. 

Two  cases  are  considered,  l.e.,  with  and  without  soil  cover  (1  ft  depth) 
over  the  floor  surface  for  radiation  protection. 

A.l  FAILURE  PROBABILITY  OF  THE  WOOD  FLOOR  SYSTEM 

This  section  presents  calculations  leading  to  the  determination  of  the 
probability  of  failure  of  the  expediently  upgraded  floor  system.  This  floor 
system  consists  of  joists,  girder,  columns,  and  stud  walls. 

A.l .1  Material  Properties 

The  entire  floor  system.  Figures  A-l  and  A-2,  consists  of  Jack  Pine 
whose  properties,  fbr  several  loading  conditions,  are  given  In  Table  A-l. 
Specific  properties  used  In  this  analysis  are  for  1  sec  load  duration.  The 
coefficient  of  variation  associated  with  each  of  these  values  Is  taken  as  0.20. 

A. 1.2  Applied  Load 

Two  load  cases  are  considered  In  the  analysis,  l.e., 

(1)  A  uniformly  distributed  mean  pressure,  p  (psl),  with  a 
coefficient  of  variation,  (Ip  ■  0.20. 

(2)  A  uniformly  distributed  mean  pressure,  'p,  plus  1  ft  of 
soil  load. 
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East  Side 
34  ft  7-1/4  inch 


North  Side 


Figure  A-l.  Basement  plan. 
(Dunes  house) 
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Figure  A-2.  Joist,  girder,  and  upper  story  partition  layout. 

(Dunes  house) 
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Designates  Partition  in  Upper  Story 


The  dead  load  of  the  floor  system  Is  neglected.  The  uniformly  dis¬ 
tributed  pressure,  j>,  Is  sssumed  to  be  applied  for  1  sec. 

A.1.3  Hamber  sizes 

Joists:  1.62S  In.  by  5.625  in.  with  an  average  spacing 

of  24.12  In. 

61rder:  5.5  In.  by  6.75  In. 

Columns:  4.0  In.  by  8.0  In. 

Studwells:  Columns  2  In.  by  4  in.  with  bracing  at  mid¬ 
height  (Figure  A-3). 

A. 1.4  Assumptions 

(1)  Joists  1  through  17  (Figure  A-2)  are  identical  and 
continuous  over  the  girder  and  the  stud  walls. 

(2)  The  connections  between  the  flooring  and  the  joists 
are  not  sufficiently  strong  to  develop  composite 
action.  Therefore,  the  joists  act  Independent  of 
the  flooring. 

(3)  The  flexibility  of  the  girder  In  calculating  joist 
stresses  Is  neglected.  The  girder  Is  assumed  to 
provide  a  rigid  support  for  the  joists. 

(4)  Resistance  along  a  member  Is  perfectly  correlated, 
l.e.,  failure  of  the  member  occurs  at  the  point 
of  maximum  load  effect. 

A.1.5  Failure  Probability  of  a  Joist 

The  joist  loading,  shear  and  bending  moment  diagrams  are  shown  in 
Figure  A-4. 

Mmax  s  M  *  9759  P  1b/1n*  (P  is  *n  lb/in.2) 

V  -  V  =  891  p  lb  (p  is  In  lb/in.2) 

max 

*  Sly  *  0.20. 

A. 1.5.1  Modes  of  Failure  -  Bending 

%  ■  ^  ■  %,  Fb  S'H 


(A-1) 
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ft-10  inch 


509p 


were  eb  *  safety  factor  In  bending 

Ng1  *  correc^on  fector  on  the  flexural  formula 
■  rupture  strength  (see  Table  A-l) 
fb  *  applied  bending  stress 
S  *  section  modulus. 


»fl  Tb  ™ 


mean  safety  factor  and  the  remaining  parameters  are 
mean  values  of  those  Identified  In  equation  (A-l). 


V  '  Vngl  +  nb  + 


where  and  are  coefficients  of  variation  of 

®b*  ^gl*  ^b*  and  respect1vely. 
a.  Uncertainty  In  the  Flexure  formula 

Due  to  the  difference  between  the  Idealized  linear  elastic  formula  and 
the  real  case,  assume  ■  0.95  with  a  uniform  distribution  between  0.90 
and  1.0,  as  shown  In  Figure  A-5  below. 


Figure  A-5.  Uniform  distribution. 
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For  a  uniform  distribution,  see  Figure  A-5. 


a  73  rrr  (Ref  22) 


(A-4) 


thus*  ngl  "  75  TM  "  °*03- 

b.  Uncertainty  In  the  Section  Modulus,  S 


s  -iEi 
*  nr 


(A-5) 


where  b  *  width  of  the  joist 
h  *  depth  of  the  joist. 

If  b  and  h  are  perfectly  correlated,  then  the  coefficient  variation 
of  S  Is 


«s  ■  V4  ♦ 4!^  ♦  <*-«) 

Assume  that  6b  *  6^  *  0.05.  Also,  due  to  humidity  effects,  let 
Ab  *  Ah  a  0.05.  Thus 

°h  ‘  nb  '  V3b  +  Ab  *  °'07 

^  -^(0.07)*  +  4(0.07)Z  +  4(0. 07 ) ( 0 . 07 )  «  0.21 

Due  to  possible  existence  of  notches,  the  mean  height  of  the  joists  is 
taken  as  0.8h.  The  section  modulus,  ?,  becomes 

; .  .  5.,„  (1,.)3 

D 
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c.  Determination  of  9b  and  flQ 
From  equations  (A-l)  and  (A-3), 


6b  =  0. 95(71 00) ( 5 .484 )/9759  p  -  3.790/p 

-  [(0.03)2  +  (0.20)2  +  (0.21 )2  +  (0.20)2]%  -  0.354 
eb 


d.  Probability  of  Failure  Due  to  Bending 


P(Fb)  1-* 


tn  (3.79/p) 
0.354 


=  1  -  $ 


[1.332  -  Jlnp  1 

- OST^J 


A. 1.5.2  Modes  of  Failure  -  Shear 


6 


v 


2Ng2  Fy  A/3V 


(A— 7) 


(A-8) 


where  0y  *  safety  factor  In  shear 

Ng2  a  correction  factor  on  the  shear  formula 
Fy  *  shear  strength  (see  Table  A-l) 
fy  *  shear  stress  due  to  applied  load 
A  ■  cross-sectional  area  of  joist 

\  '  (»gZ  *  4y  *  4  *  4,%  <A-9> 

a .  Uncertainties  In  the  Shear  Formula 

.  Thus,  flg2  "  0,03 

fl  -  0.20  (estimated) 

Fv 

■  (4  *  4  ♦ 


1 


i 


< 


f  f 


o 
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Assume  that  flg2  *  flg1 


I 


ft 


ft 


t 


a 


a 


a 


a 


a 


flA  -  [(0.07)2  +  (0.07)2  +  2(0.07)(0.07)3*  »  O.H 
*  0.20  (estimated) 

b.  Determination  of  ~9V  and 

8y  *  2tfg2  Tv  J/Zl  •  2(0.95){750)(1.625)(0.8)(5.625)/(3)<891)p 

-  3.898/p 

0  «  [(0.03)2  +  (0.20)2  +  (0.14)2  +  (0.20)2]^*  «  0.317 

°v 

c.  Probability  of  Failure  Due  to  Shear 

P(f„)  ■  1  -  ♦  -1  -♦[■  (a-,0) 

A. 1.5. 3  Joist  Failure  Probability 

The  joist  can  fail  In  flexure  or  in  shear.  As  discussed  In  Section  3.3, 
if  the  two  failure  modes  are  Independent  of  each  other  then  the  failure 
probability  of  the  joist,  P(Fj)  is 

P(Fj)  »  1  -  [1  -  P(Fb)Hl  -  P(FV)]  (A-ll) 

If  the  failure  modes  are  highly  correlated,  then 

P(Fj)  -  max[P(Fb),  P(FV)]  (A-12) 

The  actual  failure  probability  for  the  joist  Is  between  these  two 
probabilities,  thus 

max[P(Fb),  P(Fy)3<  P(Fj)  >  1  -  [1  -  P(Fb)3D  -  P(Fy)3  (A-13) 


ft 


ft 
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Failure  probabilities  of  the  joist  were  computed  and  are  shown  In 
Figure  A-6.  Failure  probabilities  P( Fb)  and  P( Fy)  were  computed  using 
equations  (A-7)  and  A-10),  respectively.  The  upper  bound  failure  probability 
was  computed  using  equation  (A-ll). 

It  is  noted  (see  Figure  A-6)  that  the  two  bounds,  l.e.,  P( Fj )  and 
P(Fv)  are  fairly  close.  In  this  analysis,  P(Fj)  Is  taken  as  the  failure 
probability  for  the  joist. 

A. 1.5. 4  Failure  Probability  of  the  Joist  System 

When  all  joists  are  Identical  and  subject  to  the  same  load  distribution 
and  Intensity,  then  conditions  between  the  joists  are  perfectly  correlated. 

On  this  basis  the  failure  probability  of  the  joist  system  Is  represented 
with  the  failure  probability  of  one  joist.  Therefore,  the  upper  bound  values 
given  In  Figure  A-6  are  conservatively  considered  as  the  failure  probability 
of  the  entire  joist  system. 

A. 1.6  Failure  Probability  of  the  Girder 

As  shown  In  Figure  A-2,  the  girder  consists  of  two  parts,  l.e.,  the 
part  between  columns  1  and  3  (part  1)  and  the  part  between  columns  3  and  the 
south  wall  (part  2).  The  two  parts  are  analyzed  separately. 

A. 1,6.1  Analysis  of  Girder,  Part  1 

The  configuration  of  this  portion  of  the  girder  Is  shown  In  Figure  A-7. 
The  loading,  P,  Is  due  to  joists. 


?  *  1278jp  (see  Figure  A-4) 

where  ?  *  mean  value  of  the  joist  load,  lb 

p  ■  mean  value  of  the  uniform  pressure,  psl. 

?  *  41.323?  *  41 .323(1278)?  *  52810. 8p 


where  M  =  mean  value  of  the  maximum  bending  moment  on  the  girder 
(at  the  center  support). 


V  -  2.502?  «  2.502(1278)p  =  3197. 6p 


where  V  *  mean  value  of  the  maximum  shear  on  the  girder 
(to  the  left  of  center  support). 


*  np  ■  o.2o 


a.  Bendinc 


\  rb  *'K 

S  .  bhf.  .  6.75(5.5)Z  .  34  03  (1n_)3 

e.  -  (0.95)(7100)(34.03)/(52810.8p)  -  4.3463/p 

’  (ng'  *  °Fh  *4*4^'  °-354 


p(Fb>  ■ 1  -  *  he 


1  -  * 


1.469  -  tni 


b.  Shear 


\  *  *g2  F,  *'3V 

8V  -  2(0.95)(750){6.75)(5.5)/3(3197.6p)  -  5.515/p 

na  •  0.317  (Determined  earlier  in  connection 
8v  with  joist  analysis) 


(A— 14) 
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P(Fy)  *1-4 


1.707  -  trip 

<>■.317“* 


(A-15) 


A. 1.6. 2  Analysts  of  Girder.  Part  2 

The  confl gu rati on  of  this  portion  of  the  girder  Is  shown  In  Figure  A-8. 
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Figure  A-8.  Girder  loading.  Part  2. 


where  ff 


where  7 


7  ■  1278p  (see  Figure  A-4) 

K  •  26.3827*  26. 382(1 278p)  *  33716.2p 
mean  value  of  the  maximum  moment  on  the  girder 
7  -  2 .2897  -  2.289(1278p)  -  2925. 3p 
mean  value  of  the  maximum  shear  on  the  girder 

*b  *  Wgl  Tb 

\  «  (0.95)(7100)(34.03)/33716.2p  •  6.808/p 

fl-  *  0.354  (computed  previously) 

^b 
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( A— 1 6) 


«%>  ■ 1  -  •  [**%&*) 

®v  "  2*g2  Tv  1/37 

0v  =  2(0.95)(750)'6.75)(5.5)/3(2925.3p)  =  6.028/p 

fl  =  0.317  (computed  previously) 
ev 

p(fvi  .  i  - 1  (u’y>)  u-17) 

Failure  probabilities  for  the  girder,  parts  1  and  2,  are  shown  In 
Figures  A-9  and  A-10.  As  previously,  three  curves  are  given,  l.e.,  probability 
of  failure  due  to  bending,  P( F^) •  probability  of  failure  due  to  shear,  P(Fy), 
and  the  upper  bound  probabilities,  P( )  and  P(Fg2)»  computed  using  equation 
(A-13).  In  each  case,  the  actual  failure  probability  Is  between  the  bounds 
of  P(Fy)  and  P(Fgi),  1  *  1,2. 

A.1.7  Failure  Probability  of  Columns 
A. 1.7.1  Existing  Columns 

The  location  and  spacing  of  columns  is  shown  in  Figures  A-l  and  A-2. 

For  a  uniformly  distributed  pressure  load  over  the  floor  surface,  the  axial 
loads  on  the  five  columns  have  the  following  values: 

P1  =  1.498P 
P2  *  4.592P 
P3  =  2.544P 
P4  *  3.655P 
P5  =  3.595P 

2 

where  P  -  1278p,  p  *  the  uniformly  distributed  floor  load  In  psl  (lb/in.  ). 

The  following  formula  (Ref  23)  was  used  for  evaluating  the  failure 
probability  of  the  five  timber  columns: 


P 

J 


0.30E 

U/d)2 


(A-l 8) 


69 


where  P  ■  column  load 

A  *  cross-sectional  area  of  column 
E  »  modulus  of  elasticity 
l  -  unsupported  column  length 

d  *  least  dimension  of  the  rectangular  cross-section  of  the  column. 

Columns  used  In  this  basement  consist  of  4  -  2  In.  x  4  In.  boards  nailed 
together.  The  columns  had  the  following  properties: 

Cross-section  ■  3.5  In.  x  6.0  In. 

Length  ■  76.0  In. 

E  *  modulus  of  elasticity  *  1.35  (10)®  psl. 

The  coefficient  of  variation  of  F„  Is  obtained  as 

c 


«...  -&f 


vce>  i&fm 


(A-19) 


where  V(  )  »  variance  of  the  given  parameter.  In  this  analysis  E  and  d 
are  considered  to  be  random  variables.  The  length,  £,  Is 
assumed  to  be  a  constant. 


(A— 20) 


ftg  ■  0.20 

nd  *  0.07  [see  tl^  In  equation  (A-6)] 
%  ■  [(0.20)2  ♦  4(0. 07)2]1*  ■  0.24 


.  i 

"gc  fc 


(A— 21 ) 


where  N  ■  correction  factor  on  the  column  fomula,  estimated  at 
9  0.95  (see  "Uncertainty  In  the  Flexure  Formula"  and 

Figure  A-5) 

f£  ■  column  stress  due  to  applied  load. 


0  *  (0.3)W_  rx/[(i/l)2^] 

c  gc  • 


(A— 22) 


where  ^  *  applied  column  load,  1  =  1 ,5 
i/d  =  76/3.5  =21.71 

0  =  (0.30)(0.95)(1,350,000)(3.5)(6.0)/[21.71)ZP1]  *  17,143/^. 


“  (fiac  +  nF  +  +  °P  ^ 

ec  9C  Fc  A  M 


(A— 23) 


n  =  0.03,  (see  A-4) 
gc 

%  -  0.24 
c 

nA  =  0.14 

-  0.20 

o  =  [(0.03)2  +  (0.24)2  +  0.14)2  +  (0.20)2]1*  -  0.34 

O  _ 


Expressions  for  computing  failure  probabilities  for  the  five  columns 


are  given 

next. 

Col  umn 

0c 

P(FC)) 

1 

1914.4? 

8.9544/? 

1  -  4  ^ 

2.1921  -  An?\ 
0.34  / 

( A— 23) 

2 

5868.6p 

2.9211/? 

1  -  4  ^ 

1 .0720  -  Anp  \ 
- 0.34  / 

(A-24) 

3 

3251.2? 

5.2727/? 

1  -  4  ^ 

1.6625  -  An?  \ 

0.S  / 

(A-25) 

4 

4671  .1? 

3.6699/? 

1  -  4  ^ 

K3T^) 

( A— 26) 

5 

4594.4? 

3.7312/? 

1  -  4  / 

*1.3167  -  An?  \ 
0.34  / 

(A-27) 

Corresponding  failure  probabilities  are  shown  In  Figure  A-ll.  P(Fr) 

Is  the  column  system  failure  probability  assuming  that  the  columns  act 
Independent  of  each  other.  It  was  computed  using  the  following  expression, 
see  equation  (49) : 


P(FJ  =  1  -  n  [1  -  P(F  *)]  (A-28) 

c  1-1  C1 

The  lower  bound  Is  represented  by  P(  Fcl>* 

A. 1.7. 2  Studwall  Columns 

The  studwall  columns  used  In  the  expedient  upgrading  of  the  basement 
are  shown  In  Figure  A-3.  See  Figure  A-l  for  the  center  to  center  dimensions. 

Studwall  column  cross-section  =  1.5  in.  x  3.5  In. 

Height,  £  =  76  In. 

£/d  -  76/(2  x  1.5)  *  25.33 

Column  load:  East  side,  Pg  *  1683p  (see  Figure  A-4) 

West  side,  Pw  *  1490p 

Studwall  columns  at  each  side  (east  or  west)  are  Identical.  Failure 
of  the  column  system  at  each  side,  therefore.  Is  represented  by  the  failure 
of  one  column. 


ec  -  TTgc  Tc  m 


1 


where  F  Is  defined  by  equation  (A-18),  thus 

C 

®c  "  °*3Ngc  E  A/CU/J)2^].  ne  «  0.34 
3  c 


For  the  east  side: 

\  •  0.3(0.95)0, 350,000)(1.5)(3.5)/[(25.33)21683p]  -  1.8706/p 
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P(F$e)  »  1  -  ♦  (A-24) 

For  the  west  side: 

6„  -  0.3(0.95)(l,350,000)(1.5)(3.5)/[(25.33)21490p]  -  2.1129/p 

v. 

p(fsw)  -i  -♦(°-7<8a.ji‘"?)  (A-25) 

Results  are  shown  In  Figure  A-12.  For  the  column  system  comprising  the 
two  studwalls,  the  lower  bound  failure  probability  Is  given  by  P(FSW).  The 
upper  bound  Is  given  by  P(F  ),  which  was  computed  using  the  expression 

P(FS)  -  1  -  [1  -  P(FSW)][1  -  P(F$e)]  (A— 26) 


A. 1.8  Failure  Probability  of  the  System 

Upper  bound  values  are  obtained  based  on  the  assumption  that  the  condi¬ 
tions  between  different  components  are  statistically  Independent.  Lower 
bound  values  are  based  on  the  perfect  correlation  assumption  between  the 
components . 

P(F*)  *  upper  bound  failure  probability 

P(F*)  -  1  -  [1  -  P(Fj)][l  -  P(Fgl )][1  -  P(Fg2)][l  -  P(FC)][1  -  P(FSW)][1  -  P(F$e)]  (A-27) 
P(F')  «  lower  bound  failure  probability 


P(F')  -  maxCP(Fj).  P(Fgl),  P(Fg2),  P(FC),  P(F$W),  P(F$e)]  (A-28) 


A. 2  PROBABILITY  OF  PEOPLE  SURVIVAL 


As  shown  In  Figure  A-13,  the  relevant  range  of  overpressures  Is  between 
0  and  about  5  psl.  In  this  range,  the  dominant  casualty  mechanism  for 
people  In  basements  Is  debris  from  the  collapse  of  the  overhead  floor  system 
and  the  upper  story.  The  upper  story  Is  expected  to  fall  In  the  range  of 
1.5  psl  to  2.9  psl.  The  following  failure  probabilities  are  extlmated  (Ref  24) 


Probability  of 


Failure  (X) 


10 

50 

90 


Probability  of  people  survival  against  structural  collapse  Is  determined 
using  the  theorem  of  total  probabilities  (Ref  25)  as 

P(SsC)  =  P(S(F)P(F)  +  P(S | F)P(F)  (A— 29) 

where  P(S  )  -  probability  of  people  survival  against  structural  collapse 

P(S|7)  *  probability  of  people  survival  given  that  the  structure 
(floor  system)  does  not  fall 

P(T)  =  probability  of  structure  survival 

P(S | F)  *  probability  of  people  survival  given  that  the  structure 
collapses  (falls) 

P(F)  =  probability  of  structural  failure  *  1  -  P(7). 

For  this  structure,  P(F)  Is  given  by  equations  (A-27)  and  (A— 28) . 

No  fatality  level  casualties  are  expected  prior  to  the  collapse  of  the 
floor  system  and,  therefore,  P(S|7)  Is  set  equal  to  1.  Probability  of 
survival  given  that  the  structure  collapses,  P(S|F),  Is  estimated  to  be  0.5. 

It  Is  based  on  the  following  reasoning. 

When  the  floor  system  over  the  basement  collapses,  the  debris  Is  not 
expted  to  affect  the  entire  shelter  area.  Several  portions  of  the  basement 
are  expected  to  be  free  of  debris.  People  located  In  these  areas  will  be 
survivors.  At  least  one  half  of  the  total  basement  area  Is  expected  to  be 
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free  of  debris  effects.  For  people  uniformly  distributed,  the  probability 
of  survival  Is,  therefore,  estimated  as  0.5. 

Probability  of  people  survival  results  are  given  In  Figure  A-14.  Two 
cases  are  considered,  i.e.,  with  and  without  soil  cover  for  radiation  pro¬ 
tection. 

The  analysis  and  results  given  here  represent  an  update  and  revision 
of  results  given  In  Reference  20. 


Overpressure,  psi 

Figure  A-14.  Probability  of  people  survival, 
upper  and  lower  bound. 


80 


APPENDIX  B 


PROBABILITY  OF  PEOPLE  SURVIVAL  IN  UPGRADED 
BASEMENTS  OF  SINGLE  FAMILY  RESIDENCES 


B.l  INTRODUCTION 

This  appendix  contains  results  on  the  probability  of  survival  of  people 
located  In  basements  of  framed,  single-family  residences  when  subjected  to 
the  blast  effects  of  a  single,  megaton  yield  unclear  weapon  detonated  near 
the  ground  surface.  Four  buildings  are  considered.  In  each,  the  basement 
Is  expediently  upgraded  against  blast  effects  by  providing  additional 
supports  for  the  joist  floor  system.  Additional  supports  are  either  stud- 
walls  or  post  and  beam  (girder  and  column)  systems.  Six  cases  are  considered: 


Building  Name 

1 .  Dunes  House 

2.  West  House 

3.  Park  House 

4.  Tea  Pot  House 


Type  of  Upgrading 

Studwall  (see  Figure  A-3) 

Girder  and  Column  (see  Figure  B-l) 

Studwall 

Studwall 

Studwall 

Girder  and  Column 


The  analysis  considering  the  Dunes  House  with  the  studwall  upgrading 
Is  described  In  Appendix  A,  which  also  contains  the  probability  of  structural 
f  llure  and  people  survival  results.  The  remaining  cases  outlined  above  are 
summarized  In  the  following  sections. 


B.2  GENERAL  ASSUMPTIONS 

(1)  The  basement  framing  system  (joists,  girders,  columns) 
for  each  building  Is  assumed  to  consist  of  "Jack  Pine" 
whose  properties  are  given  In  Table  A-l. 

(2)  The  expedient  upgrading  system  (studwall,  girder  and 
column)  Is  also  assumed  to  consist  of  Jack  Pine. 

(3)  The  upper  story  In  each  case  Is  assumed  to  fall  and  be 
removed  by  the  blast  In  the  overpressure  range  of  1.5 
to  2.9  psl  (Ref  24). 

(4)  There  Is  no  Interaction  between  the  upper  story  and  the 
basement  framing  systems,  l.e.,  the  upper  story  is  assumed 
to  cause  no  damage  to  the  basement  while  being  broken  and 
removed  by  the  blast  loading. 


Figure  B-l.  Post  and  beam  expedient  upgrading  concept. 


(5)  Dead  load  of  the  floor  system  over  the  basement  is  neglected. 
This  amounts  to  approximately  15  psf. 

(6)  People  are  assumed  to  be  uniformly  distributed  in  all 
basement  areas. 

(7)  The  only  casualty  mechanism  considered  in  the  analysis  Is 
debris  from  the  breakup  and  collapse  of  the  floor  system 
Into  the  basement  ar^a, 

(8)  Basement  walls  are  assumed  to  be  stronger  then  all  other 
structural  components  and  are,  therefore,  assumed  not  to 
fall.  Analyses  to  determine  failure  overpressures  for  the 
peripheral  basement  walls  were  net  performed.  However, 
based  on  the  results  of  full-scale  field  tests  (Ref  26) 
this  Is  a  reasonable  assumption  In  this  case. 
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8.3  OUNES  HOUSE 

The  basement  of  this  house  Is  described  In  Appendix  A  where  It  was 
analyzed  with  an  expedient  blast  upgrading  consisting  of  a  studwall  In  each 
of  the  two  joist  spans  (see  Figure  A -2).  In  this  section  the  expedient  up¬ 
grading  consists  of  the  "girder  and  column"  concept  (see  Figure  B-l)  with 
and  without  soil  cover  for  nuclear  radiation  protection.  The  size  (cross- 
section  and  length)  of  girders  used  Is  assumed  to  be  the  same  as  that  of  the 
existing  girder.  Columns  are  also  assumed  to  be  of  the  same  size  and  number 
as  the  existing  columns  and  are  assumed  to  be  Identically  spaced  and  supported. 
The  analysis  was  performed  along  the  lines  described  in  Appendix  A.  Results 
are  summarized. 

B .3.1  Failure  Probabilities 

Failure  probabilities  for  all  structural  components  making  up  the  floor 
system  and  Its  supporting  elements,  except  the  basement  walls,  are  given  In 
Figures  B-2  through  B-4.  Figure  B-2  Illustrates  failure  probabilities  for 
the  joists  and  the  existing  girder.  Failure  probabilities  for  existing 
columns  and  columns  used  with  the  expedient  upgrading  are  given  in  Figure  B-3. 
Failure  probabilities  for  the  two  sets  of  girders  used  in  the  expedient  up¬ 
grading  are  given  in  Figure  B-4.  Each  of  the  curves  is  an  upper  bound  on  the 
particular  failure  probability  and  was  determined  using  equation  (49).  The 
bounds  on  the  probability  of  failure  of  this  expediently  upgraded  floor  system 
are  given  In  Figure  B-5  for  the  case  without  soil  cover.  In  this  figure, 

P(F*),  the  upper  bound.  Is  based  on  equation  (49)  and  P(F'),  the  lower  bound, 
represents  the  failure  probability  of  upgrading  column  2  located  In  the  east 
joist  span  (see  Figure  A-2). 

B.3.2  People  Survival  Probabilities 

Probabilities  of  people  survival  are  given  in  Figure  B-6,  which  Includes 
two  cases,  l.e.,  with  and  without  sold  cover.  Probability  of  survival  Is 
against  the  effects  of  debris  produced  by  the  breakup  of  the  floor  system. 

B.4  WEST  HOUSE 

This  an  existing  single-family  dwelling  whose  basement  floor  plan  Is 
shown  in  Figure  B-7.  The  floor  system  over  the  basement  consists  of  a  subfloor 
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existinq  qirder  failure  Probabilities,  Dunes  House. 


failure  probabilities.  Dunes  House. 
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Without  Soil  Cover 


Probability  of  people  survival,  upper  and  lower  bound.  Dunes  House 


in.  Block  Mall 


West  House,  plan. 


and  a  finish  floor  supported  by  2  In.  x  10  In.  joists  spaced  at  16  In.  The 
two  joist  spans  are  simply  supported.  The  existing  girder  Is  6  In.  x  10  In. 
and  consists  of  two  parts.  Part  1  extends  from  the  front  wall  of  the  house 
to  column  2.  Part  2  extends  from  column  2  to  the  rear  wall  of  the  house. 

The  three  columns  have  a  cross-section  of  6  In.  by  6  In.  Their  unsupported 
length  Is  77  In. 

The  floor  system  Is  assumed  to  be  upgraded  using  studwalls  In  each  of 
the  two  joist  spans.  The  studs  are  2  In.  x  4  In.  and  are  spaced  at  16  In. 
as  are  the  joists.  Their  total  height  Is  70  In.  They  are  braced  at  half¬ 
height  as  shown  In  Figure  A-3.  The  analysis  was  performed  along  the  lines 
described  In  Appendix  A.  Results  are  summarized. 

B .4.1  Failure  Probabilities 

Failure  probabilities  for  the  joists  and  the  two  girders  are  given  in 
Figure  B-8.  Failure  probabilities  for  the  columns  and  studwalls  are  given 
In  Figure  B-9.  Each  of  these  curves  is  an  upper  bound  on  the  particular 
probability  and  was  determined  using  equation  (49). 

The  bounds  on  the  probability  of  failure  of  the  whole  floor  system. 
Including  columns  and  studwalls,  are  given  In  Figure  B-10.  P(F)*,  the  upper 

bound  probability  of  failure  of  the  system,  was  obtained  using  equation  (49) 
and  P(F)',  the  lower  bound,  represents  the  failure  probability  of  the  stud¬ 
walls. 

B.4.2  People  Survival  Probabilities 

Probabilities  of  people  survival  against  the  effects  of  debris  from  the 
collapse  of  the  floor  system  Into  the  basement  are  given  In  Figure  B-ll.  Two 
cases  are  considered,  l.e.,  with  and  without  soil  cover  for  radiation  protec¬ 
tion.  Probability  of  survival  Is  against  the  effects  of  debris  produced  by 
the  collapse  of  the  floor  system  Into  the  basement. 

B.5  PARK  HOUSE 


This  Is  an  existing  residence  whose  basement  floor  plan  Is  shown  In 
Figure  B-12.  The  floor  system  over  the  basement  consists  of  a  subfloor  and 
a  finish  floor  supported  by  2  In.  x  8  In.  joists  spaced  at  12  In.  on  center 
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Probability  of  floor  system  failure,  upper  and  lower  bound,  West  House 
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Figure  B-12.  Park  House,  plan. 


The  joists  are  continuous  over  the  centrally  located  girder  and  columns.  The 
girder  Is  8  In.  x  9  In.  and  consists  of  three  separate  parts.  The  first  part 
spans  from  the  east  wall  to  column  2  (see  Flg-rz  B-12),  the  second  part  spans 
from  column  2  to  column  4,  and  the  third  part  spans  from  column  4  to  the  west 
wall.  The  five  wood  columns  have  an  cross-section  of  6  In.  x  6  In.  Their 
unsupported  length  Is  approximately  7  ft  83s  In. 

The  floor  system  Is  assumed  to  be  expediently  upgraded  using  studwalls 
In  each  of  the  two  joist  spans.  This  expedient  upgrading  Is  Illustrated  In 
Figures  B-13  and  B-14.  The  studs  are  2  In.  x  4  In.  and  are  spaced  at  12  In., 
I.e.,  one  under  each  joist.  Their  total  height  Is  shown  In  Figure  B-14  and 
they  are  braced  at  half-height.  In  addition  to  structurally  upgrading  the 
floor,  the  basement  shelter  Is  also  assumed  to  be  mounded  with  soil  up  to  the 
top  of  the  floor  as  shown  in  Figure  B-13. 

The  analysis  was  performed  along  the  lines  described  In  Appendix  A. 
Results  are  summarized. 

B .5.1  Failure  Probabilities 

Failure  probabilities  for  the  joists  and  the  three  timber  girders 
supporting  them  are  shown  in  Figure  B-15.  The  combined  failure  probabilities 
for  the  columns  are  given  in  Figure  B-16  together  with  the  failure  probability 
of  the  studwalls.  Each  of  these  curves  represents  an  upper  bound  and  was 
determined  using  equation  (49). 

The  bounds  on  the  probability  of  failure  of  the  whole  floor  system. 
Including  columns  and  studwalls,  are  given  In  Figure  B-17.  P(F)*,  the  upper 

bound  probability  of  failure  of  the  system,  was  obtained  using  equation  (49). 
P(F)',  the  lower  bound  for  the  system.  Is  the  failure  probability  of  the 
studwal 1 s . 


B.5.2  People  Survival  Probabilities 

Probabilities  of  people  survival  against  the  effects  of  debris  from  the 
collapse  of  the  floor  system  into  the  basement  are  given  In  Figure  B-18. 

Two  cases  are  considered,  I.e.,  with  and  without  soil  cover  for  radiation 
protection. 


Figure  B-13.  Elevation  (Section  A-A). 
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Figure  B-14.  Studwall  expedient  upgrading. 
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Failure  probabilities  of  columns  and  studwalls.  Park  House. 


Drobability  of  people  survival,  upper  and  lower  bound,  Dark  House. 


B.6  TEA  POT  HOUSE 


This  Is  a  two-story,  single-family  residence  originally  constructed  and 
tested  In  Nevada  (Ref  26).  This  building  has  a  full  basement  with  a  back 
entrance,  an  entrance  from  the  house,  and  six  window  wells.  The  basement 
plan  Is  shown  In  Figure  B-19. 

The  floor  over  the  basement  consists  of  a  subfloor  and  a  finish  floor 
supported  by  2  In.  x  8  In.  joists  spaced  at  16  In.  centers.  The  joists 
(assumed  to  be  continuous  over  the  33  ft  4  In.  length  of  the  house)  are 
supported  by  two  6  In.  x  8  In.  girders  and  the  basement  walls.  The  two 
girders  are  supported  by  four  steel  pipe  columns  and  the  basement  walls.  The 
peripheral  basement  walls  are  made  of  concrete  block.  Two  expedient  upgrading 
schemes  were  considered  and  are  described  as  follows: 

a.  Scheme  1 

The  two  long  (13  ft  4  In.)  joist  spans  were  each  assumed  to  be  supported 
by  a  studwall  located  halfway  between  the  columns  and  the  walls.  This  con¬ 
cept  calls  for  a  2  In.  x  4  In.  stud  under  each  joist. 

Entranceways  Into  the  basement  are  assumed  to  be  closed  (blocked)  by 
means  of  expedient  blast  closures.  Window  glass  Is  assumed  to  be  removed  and 
the  openings  are  also  assumed  to  be  blocked  by  means  of  expedient  blast 
closures. 

The  basement  Is  mounded  with  soil  on  the  outside  up  to  the  first  floor 
level,  about  2  ft.  One  foot  of  soil  is  assumed  to  be  placed  on  the  first 
floor  for  fallout  radiation  protection. 

b.  Scheme  2 

This  expedient  upgrading  scheme  Is  the  same  as  the  first  scheme  except 
that  Instead  of  a  studwall,  the  two  joist  spans  are  assumed  to  be  upgraded 
by  girders  and  columns  located  halfway  between  the  existing  columns  and  the 
walls.  The  girder  Is  assumed  to  be  of  the  same  size  and  the  same  material 
as  the  existing  girder.  The  columns  consist  of  "Southern  Pine,"  have  a 
6  In.  x  4  In.  cross-section,  are  8  ft  long,  and  have  the  following  properties 
with  respect  to  an  axial  load: 


Fc  (compression  parallel  to  the  grain)  *  1350  psl 

E  (modulus  of  elasticity)  =  1,700,000  psl. 

Upgrading  columns  are  assumed  to  have  the  same  spacing  as  the  existing 
columns.  This  concept  is  Illustrated  In  general  in  Figure  B-l. 

B.6.1  Failure  Probabilities,  Scheme  1,  Studwall  Upgrading 

Failure  probabilities  for  the  joists  and  existing  girders  are  shown  In 
Figure  B-20.  Failure  probabilities  for  existing  columns  and  the  studwalls 
used  for  upgrading  are  shown  in  Figure  B-21 .  Upper  and  lower  bounds  on  the 
failure  probability  of  the  system  as  a  whole  is  shown  in  Figure  B-22.  In 
this  case  the  lower  bound  is  the  failure  probability  of  the  studwalls  also 
shown  In  Figure  B-21.  The  upper  bound  was  computed  using  equation  (49). 

B.6.2  People  Survival  Probabilities,  Scheme  1,  Studwall  Upgrading 

People  survival  probabilities  are  presented  in  Figure  B-23  and  include 
two  cases,  l.e.,  with  and  without  soil  cover  for  fallout  radiation  protection. 
Probability  of  survival  Is  against  the  effects  of  debris  produced  by  the 
breakup  of  the  floor  system  over  the  basement. 

B.6.3  Failure  Probabilities,  Scheme  2,  Girder  and  Column  Upgrading 


Failure  probabilities  for  the  joists,  the  existing  girders,  and  the 
girders  used  in  the  expedient  upgrading  are  given  in  Figure  B-24.  Failure 
probabilities  for  existing  columns  and  the  columns  used  in  the  expedient  up¬ 
grading  are  shown  In  Figure  B-25.  Upper  and  lower  bounds  on  the  failure 
probability  of  the  system  as  a  whole  are  shown  in  Figure  B-26.  In  this  case 
the  lower  bound  is  the  failure  probability  of  the  columns  used  in  the  expedient 
upgrading.  This  Is  also  shown  In  Figure  B-25.  The  upper  bound  was  computed 
using  equation  (49), 

B.6.4  People  Survival  Probabilities,  Scheme  2, 

C4rder  and  Column  Upgrading _ 

People  survival  probabilities  for  this  concept  are  given  In  Figure  B-27. 
Two  cases  are  considered,  l.e.,  with  and  without  soil  cover  for  fallout 
radiation  protection.  Probability  of  survival  Is  against  the  effects  of 
debris  produced  by  the  breakup  of  the  floor  system  over  the  basement.  It 
Is  evident  that  the  difference  between  the  upper  and  the  lower  bounds  on  the 
probability  of  survival  Is  negligible  In  this  case. 


s  and  existinq  qirders  failure  probabilities.  Scheme  1,  T:a  Pot  House 


Studwal Is  and  existing  columns  failure  probabilities.  Scheme  1  Tea  Pot  House 


Finure  B-24.  Joist  and  qirder  failure  probabilities.  Scheme  2,  Tea  Pot  House. 


Failure  probabilities  of  col  wans.  Scheme  2,  Tea  Pot  House 


Probability  of  floor  system  failure,  UDDer  and  lower  bound.  Scheme  2,  Tea  Pot  House. 


Fiqure  B-27.  Probability  of  people  survival,  upoer  and  lower  bound.  Scheme  2,  Tea  pot  House 
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APPENDIX  C 

STRUCTURAL  FAILURE  AND  PEOPLE  SURVIVAL  PROBABILITY  DATA 

This  appendix  contains  detailed  results  on  the  probability  of  structural 
failure  and  the  probability  of  people  survival  for  the  reinforced  concrete 
shelters  described  in  Chapter  4.  The  general  concept  of  the  basic  basement 
shelters  is  illustrated  in  Figure  1.  These  basements  were  designed  for  live 
loads  in  the  range  from  50  psi  to  250  psl  and  span  lengths  from  12  ft  to  20  ft. 
The  basic  design  data  are  given  in  Table  2.  Each  of  the  12  slabs  was  analyzed 
as  upgraded  using  four  expedient  upgrading  schemes  illustrated  in  Figure  9. 

This  resulted  in  60  sets  of  shelters  whose  analysis  data  are  included  in 
Table  3.  Results  included  here  are  upper  and  lower  bounds  on  the  probability 
of  slab  failure  and  upper  and  lower  bounds  on  the  probability  of  people  survival. 
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Figure  C-48.  Probability  of  people  survival  (upper  and  lower  bounds)  case  5D 
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Figure  C-120.  Probability  of  people  survival  (upper  and  lower  bounds)  case  12E 
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